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ABSTRACT 


The recent discovery by Cantalupo et al. (2014) of the largest (~ 500 kpc) and luminous 
(L ~ 1.43 X 10^^ erg s“^) Lya nebula associated with the quasar UM287 {z = 2.279) poses a great 
challenge to our current understanding of the astrophysics of the halos hosting massive z ^ 2 galaxies. 
Either an enormous reservoir of cool gas is required M ~ 10^^ exceeding the expected baryonic 
mass available, or one must invoke extreme gas clumping factors not present in high-resolution cos¬ 
mological simulations. However, observations of Lya emission alone cannot distinguish between these 
two scenarios. We have obtained the deepest ever spectroscopic integrations in the He II A1640 and 
C IV A1549 emission lines with the goal of detecting extended line emission, but detect neither line to 
a 3cr limiting SB~ 10“^®ergs“^ cm“^ arcsec“^. We construct simple models of the expected emission 
spectrum in the highly probable scenario that the nebula is powered by photoionization from the cen¬ 
tral hyper-luminous quasar. The non-detection of Hell implies that the nebular emission arises from 
a mass Me < 6.4 x 10^° Mq of cool gas on ~ 200 kpc scales, distributed in a population of remarkably 
dense (nn ^ 3 cm“^) and compact {R < 20 pc) clouds, which would clearly be unresolved by current 
cosmological simulations. Given the large gas motions suggested by the Lya line {v ~ 500km/s), it 
is unclear how these clouds survive without being disrupted by hydrodynamic instabilities. Our work 
serves as a benchmark for future deep integrations with current and planned wide-held IFU spec¬ 
trographs such as MUSE, KGWI, and KMOS. Our observations and models suggest that a ~ 10 hr 
exposure would likely detect ~ 10 rest-frame UV/optical emission lines, opening up the possibility of 
conducting detailed photoionization modeling to infer the physical state of gas in the circumgalactic 
medium. 

Subject headings: galaxies: formation — galaxies: high-redshift — intergalactic medium — circum¬ 
galactic medium 


1. INTRODUCTION 

In the modern astrophysical lexicon, the intergalactic 
medium (IGM) is the diffuse medium tracing the large- 
scale structure in the Universe, while the so-called cir¬ 
cumgalactic medium (CGM) is the material on smaller 
scales within galactic halos (r < 200 kpc), for which non¬ 
linear processes and the complex interplay between all 
mechanisms that lead to galaxy formation take place. 


Whether one is studying the IGM or the GGM, for 
decades the preferred technique for characterizing such 
gas has been the analysis of absorption features along 


background sightlines (e.g. 

Croft et al. 2002 Bergeron 

et al. 

120041 IHennawi et al. 


Hennawi & Prochaska 

2007f 

Prochaska & Hennawi|2009t 

Rudie et al.|2012t|Hen- 
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nawi fc Prochaska1|2013[ [Farina et al.|[20T^ [Prochaska] 
et al.||2013a|b| |Lee et al.||2014p. However, as the ab¬ 
sorption studies are limited by the rarity of suitably 
bright background sources near galaxies, and to the one¬ 
dimensional information that they provide, they need 
to be complemented by the direct observation of the 
medium in emission. 

In particular, it has been shown that UV background 
radiation could be reprocessed by thes e media and be 
detectable as fluorescent Lya emission (|Hogan fc 


mann|I987l Binette et al.|I993[ Gould fc Weinberg|1996[ 

Cantalupo et al. 2005). However, current facilities are 

still not capable of revealing such low radiation levels, 
e.g. an expected surface brightness (SB) of the o rder of 
SB t,v» ~ 1 0~^° ergs“^ cm“^ arcsec”^ (see e.g. Rauch 


et al. 


2008). Nonetheless, this signal can be boosted to 


observable levels by the intense ionizing flux of a nearby 


surroundings 1 

Rees|1988 iHaiman & Rees 12001 Alam & 

Miralda-Escud 

e 2002 Cantalupo et al. MIS), shedding 


Detecting this fluorescence signal has been a subject 
of significant interest, and several studies which specif- 


imity to a quasar (e.g., Fynbo et al. 1999 

Francis & 

Bland-Hawt 

lorn 20041 Cantalupo et al.||2007[ |Rauch 

et al. 2008| 

Hennawi & Prochaska 2013p nave so far a 

not straight 

forward interpretation. But recently |Can- 

talupo et al 

(2012) identified a population of compact 
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Lya emitters with rest-frame equivalent widths exceed¬ 
ing the maximum val ue expected from star -formation, 
> 240A (e.g., Chariot & Fall 1993), which are 


the best candidates to date tor fluorescent emission pow¬ 
ered by a proximate quasar. 

Besides illuminating nearby clouds in the IGM, a 
quasar may irradiate gas in its own host galaxy or CGM. 
A number of studies have reported the detection of ex- 
tende d Lya emission in the vicinity oi z 2 — 4 quasars 
(e .g., Hu & Cowiep987 Heckman et al ||1991a|b| Chris- 


terisen et al.||2UU'6[ |Hermawi et al.||2l)l)9| [hjorth et al. 

2012p, but detailed comparison is hampered by the differ¬ 


ent methodologies of these studies. Although extended 
Lya nebulae on scales of ~ 100 kpc (up to 250 kpc) 


(HzRGs; e.g. McGarthy||1993 van Ojik et al.|1997l iReu- 

land et al.l 

20031 IVillar-Martin et al.||2003b|al IReuland 

0 

0 

Miley & De Breuck|2008 ), these objects have 


the powerful radio jets and the ambient medium, com¬ 
plicating the interpretation of the observations. Nebulae 
of comparable size and luminosity have similarly been ob- 


blobs’ (LABs) (e.g., 

Steidel et al. 2000[ 

Matsuda et al. 

2004 

Uey et al. |2005t|bmith &: Jarvis||20( 

J71 IGeach et al. 

2007 

Prescott et al. 1120091 1 Yang et al. 

120111120121120141 

Arrigoni IJattaia et ai.||2014bp which o 

0 not show direct 


the L ABs are also frequently associated with obscured 


AGN ( Geach et al. 2009 Overzier et al. 2013 Prescott 


let al.l 


_ (although lacking powerful radio jets), the 

mechanism powering their emission remains controver¬ 
sial with at least four proposed which may even act to- 


(Geach et al. 2009 

Overzier et al. 2013), (ii) shock heated 

gas by galactic su] 
cooling radiation f 

^erwinds ('I'aniguchi et al. 

2001), (hi) 

rom cold-mode accretion (e.g., h'ardal 

et al. 120011 lYang et al. 120061 IFaucher-Giguere et a .pOlOl 

Rosdahl & Blaizot||2012p, and (iv) resonant scattering of 


Hayes et al. ][ml| |Gen fc ZhengllMdl r 


i'he largest and most luminous Lya nebula known is 
that around the quasar UM 287 (Amag=17.28) at z = 
2.279, recently discovered by Gantalupo et al.|(|2014) in 


a narrow-band imaging survey o t hyper-luminous quasars 


(Arrigoni Battaia et al. |2014a|. Its size of 460 kpc and 


average Lya surface brightness of SBi^ya = 6.0 x 10“^® 
ergs“^ cm“^ arcsec“^ (from the 2cr isophote), which cor¬ 
responds to a total luminosity oiL^ya — (2.2±0.2) x 10'^'^ 
erg s“^, make it the largest reservoir (Me ~ 10^^ of 
cool (T ~ 10^ K) gas ever observed around a QSO. The 
emission has been explained as recombination and/or 
scattering emission from the central quasar and has been 
regarde d as the hrst direct det ection of a cosmic web fil¬ 
ament (Gantalupo et al.||2014|). 


However, as discussed m Gantalupo et al. (2014), Lya 
emission alone does not allow to break the degeneracy 
between the dumpiness or density of the gas, and the 
total gas mass. Indeed, in the scenario where the nebula 
is ionized by the quasar radiation, the total cool gas mass 
scales as Me ^ 10^^(7“^/^ Mq, where C = (n%,)l(n^Y is 
a clumping factor introduced by Gantalupo et al. (2014) 
to account for the possibility of higher density gas unre¬ 


solved by the cosmological simulation used to model the 
emission. Thus, if one assume (7=1, the implied cool 
gas mass in the extended nebula is exceptionally high for 
the expected dark matter halo inhabited by a z 'y 2 — 3 
quasar, i.e. Mjjm = 10^^ ® Mq (|White et al.|2012|). This 
is further aggravated by the fact that current cosmolog¬ 
ical simulations show that only a small fraction (^ 15%, 


Fumagalli et al.||20141 |Faucher-Giguere et al.||201^ |Can-| 

talupo et al.||2014 ) of the total baryons reside in a phase 

(T < 5 X 10"^ K) sufficiently cool to emit in the Lya line. 
A possible solution to this discrepancy would be to as¬ 
sume a very high clumping factor up to (7 ~ 1000, which 
would then imply a large population of cool, dense clouds 
in the CGM and extending into the IGM, wh ich are un¬ 
resolved by current cosmological simulations (Gantalupo 

eFalllMdl ). 


Both ot these scenarios, whether it be too much cool 
gas, or a large population of dense clumps, are reminis¬ 
cent of a similar problem that has e merged from absorp¬ 
tion line studies of the quasar CGM (Hennawi et al. 2006 


Hennawi &: Prochaska|2007l [Frochaska fc Hennawi 2009. 

Prochaska et al. 2013a|b 2014p . This work reveals sub- 

stantial reservoirs of cool gas > 10^° Mq, manifest as a 
high covering factor ~ 50% of optically thick absorption, 
several time s larger than predicted by hydrodynamical 
simul ations (|Fumagalli et al.|2014||Faucher-Giguere et ahj 
2014). This conflict most likely indicates that current 
simulations fail to capture essential as pects of the hy¬ 
drodynamics in massive halos at x ~ 2 (Prochaska et al.| 
2013b Fumagalli et al.||2014), perhaps tailing to resolve 
the formation of clumpy structure in cool gas, which in 
the most extreme cases give rise to giant nebulae like 
UM287. 

In an effort to better understand the mechanism pow¬ 
ering the emission in UM287, and further constrain 
the physical properties of the emitting gas, this paper 
presents the result of a sensitive search for emission in 
two additional diagnostics, namely He II A1640A and 
G IV AI 549 A. The detection of either of these high- 
ionization emission lines in the extended nebula, would 
indicate that the nebula is ‘illuminated’ by an intense 
source of hard ionizing photons E > 4Ryd, and would 
thus establish that photoionization by the quasar is the 
primary mechanism powering the giant Lya nebula. As 
we will show in this work, in a photoionization scenario 
where He II emission results from recombinations, the 
strength of this line is sensitive to the density of the gas in 
the nebula, which can thus break the degeneracy between 
gas density and gas mass described above. In addition, 
because He II is not a resonant line, a comparison of its 
morphology and kinematics to the Lya line can be used 
to test whether Lya photons are resonantly scattered. 
On the other hand, a detection of extended emission in 
the G IV line can provide us information on the metal- 
licity of the gas in the GGM, and simultaneously con¬ 
strain the size at which the halo is metal-enriched. To 
interpret our observational results, w e explo it the models 
presented by Hennawi & Prochaska (2013) and already 
us ed in the context ot extended Ly a nebulae (i.e. LABs) 
in Arrigoni Battaia et al. (2014b), and show how a sen- 


^ The He II A1640A is the first line of the Balmer series emitted 
by the Hydrogen-like atom He+, i.e. corresponding to the Ho line. 
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sitive search for diffuse emission in Lya and additional 
diagnostic lines can be used to constrain the physical 
properties of the quasar GGM. 

This paper is organized as follows. In ^ we describe 
our deep spectroscopic observations, the Sta reduction 
procedures, and the surface brightness limits of our data. 
In we present our constraints on emission in the C IV 
anoHe II lines, and our analysis of the kinematics of the 
Lya line. In ® we present the photoionization models 
for UM287 andin ^we compare them with our observa¬ 
tional res ults a nd with absorption spectroscopic studies 
( |5.I| and 15.2). In ^we discuss which other lines might 
be observable with current facilities. In 0 we further 
discuss some of the assumptions made in our modeling. 
Finally, ^summarizes our conclusions. 

Throu^out this paper, we adopt the cosmological pa¬ 
rameters Hq = 70 km s“^ Mpc“^, Qm = 0.3 and 
Da = 0.7. In this cosmology, I" corresponds to 8.2 phys¬ 
ical kpc at z = 2 .279. All magnitudes are in the AB 
system ( Oke||1974 ). 


2. OBSERVATIONS AND DATA REDUCTION 
Two moderate resolution (FWHM~ 300 km s“^) spec- 


Oke et al. 

1995) 

g 4, in mu 

iti-slit 


mode with custom-designed slitmasks. We used the 600 
lines mm“^ grism blazed at 4000 A on the blue side, re¬ 
sulting in wavelength coverage of « 3300 — 5880 A, which 
allows us to cover the location of the C IV and He II lines. 
The dispersion of this grism is ~ 4 A per pixel and our 
I"slit give a resolution of FWHM ~ 300 km s“^. We 
observed each mask for a total of ~ 2 hours in a series of 
4 exposures. 

Figure H shows the position of the two I"-slits (red 
and bluejon top of the narrow-band image (matching 


Cantalupo et al. 

(2014D 

Gantalupo et al. 

(2014 

21.54 ± (J.(J6) radi 

o-loud 


. We remind the reader that 
found a optically faint (V = 


redshift, and at a projected distant of 24.3 arcsec (^ 200 
kpc) from the bright UM287 quasar (‘QSO a’). The first 
slit orientations was chosen to simultaneously cover the 
extended Lya emission and the UM287 quasar (blue slit), 
whereas the second (red slit) was chosen to cover the 
companion quasar ‘b’ together with the diffuse nebula. 
By covering one of the quasars with each slit orienta¬ 
tion we are thus able to cleanly subtract the PSF of the 
quasars from our data (see Section 

The 2-d spectrosco pic data reduction is pe r forme d ex¬ 
actly as described in Hennawi & Prochaska (2013) and 
we refer the reader to that work for additional details. In 
what follows, we briefly summarize the key elements of 
the data reduction procedure. All data were reduced us¬ 
ing the LowRedux pipelin^ which is a publicly available 
collection of custom codes written in the Interactive Data 
Language (IDL) for reducing slit spectroscopy. Individ¬ 
ual exposures are processed using standard techniques, 
namely they are overscan and bias subtracted and flat 
fielded. Cosmic rays and bad pixels are identified and 
masked in multiple steps. Wavelength solutions are de¬ 
termined from low order polynomial fits to arc lamp spec- 

® http:// WWW. ucolick.org/~xavier/LowRedux 


tra, and then a wavelength map is obtained by tracing 
the spatial trajectory of arc lines across each slit. 

We then perform the sky and PSF subtraction as a 
coupled problem, using a no vel custom algorithm that we 
briefly summarize here (see Hennawi & Prochaska 2013 
for additional details). We adopt an iterative procedure, 
which allows us to obtain the sky background, the 2-d 
spectrum of each object, and the noise, as follows. Fir^ 
we identify objects in an initial sky-subtracted imag^ 
and trace their trajectory across the detector. We then 
extract a 1-d spectrum, normalize these sky-subtracted 
images by the total extracted flux, and fit a B-spline pro¬ 
file to the normalized spatial light profile of each object 
relative to the position of its trace. Given this set of 2-d 
basis functions, i.e. the flat sky and the object model 
profiles, we then minimize chi-squared for the best set 
of spectral B-spline coefficients which are the spectral 
amplitudes of each basis component of the 2-d model. 
The result of this procedure are then full 2-d models 
of the sky-background, all object spectra, and the noise 
(tr^). We then use this model sky to update the sky- 
subtraction, the individual object profiles are re-fit and 
the basis functions updated, and chi-square fitting is re¬ 
peated. We iterate this procedure of object profile fitting 
and subsequent chi-squared modeling four times until we 
arrived at our final models. 

For each slit, each exposure is modeled according to 
the above procedure, allowing us to subtract both the 
sky and the PSF of the quasars. These images are regis¬ 
tered to a common frame by applying integer pixel shifts 
(to avoid correlating errors), and are then combined to 
form final 2-d stacked sky-subtracted and sky-and-PSF- 
subtracted images. The individual 2-d frames are op¬ 
timally weighted by the (S/N)^ of their extracted 1-d 
spectra. The final result of our data analysis are three 
images: 1) an optimally weighted average sky-subtracted 
image, 2) an optimally weighted average sky-and-PSF- 
subtracted image, and 3) the noise model for these im¬ 
ages cr^. The final noise map is propagated from the indi¬ 
vidual noise model images taking into account weighting 
and pixel masking entirely self-consistently. 

Finally, we flux calibrate our d a ta fol lowing the pro¬ 
cedure in Hennawi & Prochaska (|2013 ). As standard 
star spectra were not typically taken immediately be¬ 
fore/after our observations, we apply an archived sen¬ 
sitivity function for the LRIS B600/4000 grism to the 
1 -d extracted quasar spectrum for each slit, and then in¬ 
tegrate the flux-calibrated spectrum against the SDSS 
g-band filter curve. The sensitivity function is then 
rescaled to yield the correct SDSS g-band photometry. 
Since the faint quasar is not clearly detected in SDSS, 
we only used the g-band magnitude of the UM287 quasar 
to calculate this correction. Note that this procedure is 
effective for point source flux-calibration because it al¬ 
lows us to account for the typical slit-losses that affect a 
point source. However, this procedure will tend to un¬ 
derestimate our sensitivity to extended emission, which 
is not affected by these slit-losses. Hence, our procedure 
is to apply the rescaled sensitivity functions (based on 
point source photometry) to our 2-d images, but reduce 
them by a geometric slit-loss factor so that we properly 

® By construction, the sky-background has a flat spatial profile 
because our slits are flattened by the slit illumination function. 
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iri 



RA (J2000) 

Figure 1. 10-hours narrow-band image matchin g the Lva line at 
the re dshift of UM287 (adapted from Figure 1 of|Cantalupo et al.| 
|2014||. ‘QSO a’ is the quasar UM287, while ‘QSU b' is the taint 
companion quasar . The red and blue lines highlight the position 
of the l"slits chosen to study the extended emission in this work. 
Note that a Lyman Alpha Emitter (‘LAE’) and a continuum source 
(‘C’) fall within the ‘red’ slit (see Figure]^. 


treat extended emission. To compute the slit-losses we 
use the measured spatial FWHM to determine the frac¬ 
tion of light going through our 1.0" slits, but we do not 
model centering er rors (see Section for a tes t of our 
calibration, and see |Hennawi fc Prochaska|2013 for more 
details). 

Given this flux calibration, the la SB limit of our ob¬ 
servations are SBi^- = 1.3 x 10“^® ergs“^ cm“^ arcsec”^, 
and SBio- = 1.5 x 10“^® ergs“^ cm“^ arcsec”^ for C IV 
and He II, respectively. This limits are obtained by aver¬ 
aging over a 3000 km s“^ velocity interval, i.e. ±1500 km 
s“ ^ on either side of the systemi c redshift of the UM28 7 
quasar, i.e. z = 2.279 ± 0.001 (|McIntosh et al.||1999D, 
at the C IV and He II locations, and a 1" x I" aper- 
tura^ This limits (approximately independent of wave- 
len gth) are about 3x the Icr limit in 1 arcsec^ quoted 
by Cantalupo et al.| (20I4) for their ~ 10 hours narrow- 
band exposure targeting the Lya line, i.e. 5 x 10“^® 
ergs“^ cm“^ arcsec“^ Note that we choose this ve¬ 
locity range to enclose ml the extended Lya emission, 
even after smoothing (see next Section j|^ and b ecause 
the narrow-band image of Cantalupo et m.|(2014) covers 


approximately this width, i.e. Av ^ 2400 km s“^. 

Further, it is important to stress here that the line 
ratios we use in this work are only from the spectroscopic 
data (we do not use the NB data for the Lya line), and 


Obviously, if we use a smaller velocity aperture we get a more 
sensitive limit, i.e. SSumit = ^ e.g. we ob¬ 

tain SBia = 7.3 x 10“^® ergs“^ cm“^ arcsec“^ for a 700 km s“^ 
velocity interval. 

Note that spatial averaging allow us to achieve more sensitive 
limits. If we consider an aperture of l"x20”, we reach = 

3.7 x 10“^® ergs“^ cm“^ arcsec“^ at the location of the Lya line. 


hence they are independent of any errors in the absolute 
calibration. Although we do not use the NB data in our 
analysis, we show in the next section that our results 
are consistent with the NB imaging, and thus robustly 
calibrated. 


3. OBSERVATIONAL RESULTS 


Following Hennawi & Prochaska (2013), we search for 
extended Lya, C IV, and He II emission by constructing 
a X image 




= E 


(DATAi - MODEL,) 


( 1 ) 


where the sum is taken over all Alpix pixels in the image, 
‘DATA’ is the image, ‘MODEL’ is a linear combination of 
2-d basis functions multiplied by B-spline spectral ampli¬ 
tudes, and cr is a model of the noise in the spectrum, i.e. 
0-2 = SKY ± OBJECTS TREADNOISE. The ‘MODEL’ 
and the cr^ are obtained du ring our data reduction pro- 
cedure (see Section © and |Hennawi & Prochaska||2013| 
for details). 

Eigures and show the two-dimensional spectra for 
the slits in Eigurej^ plotted as y-maps. Note that if our 
noise model is an accurate description of the data, the 
distribution of pixel values in the X'Hraps should be a 
Gaussian with unit variance. In these images, emission 
will be manifest as residual flux, inconsistent with be¬ 
ing Gaussian distributed noise. The bottom row of each 
figure shows the Xsky map (only sky subtracted) at the 
location of the Lya, C IV, and He II, respectively. Even 
in these unsmoothed data the extended Lya emission is 
clearly visible up to ^ 200 kpc (^ 24") from ‘QSO b’, 
along the ‘red’ slit (Eigurel^. This emission has SEpya = 
(6.3 ± 0.4) X 10“^® ergs“^cm“^ arcsec”^, calculated in 
a I" X20"apertur and over a 3000 km s“^ velocity in¬ 
terval (blue box in Figure 2). This value is in agreement 
with the emission detected in the continuum-subtracted 


image presented in Cantalupo et al. (2014) within a 
l"x20"aperture at the same position within the slit, 
i.e. SBpyQ = (7.0 ± 0.1) x 10“^® ergs“^ cm“^ arcsec”^. 
Along the ’blue’ slit, the extended emission is inevitably 
mixed with the PSF of the hyper-luminous UM287 QSO, 
making PSF subtraction much more challenging. Never¬ 
theless, we compute the emission in the extended Lya 
line in an aperture of about l"xl3"aperture (from 40 
to 150 kpc) and within 3000 km s“^, after subtract¬ 
ing the PSF of the quasar (see Figure J^. Again, 
we find that surface brightness measuredfrom spec¬ 
troscopy (SBpyct = 1.4 X 10“^”^ ergs“^ cm“^ arcsec”^), 
and from narrow-band imaging (SBpya = 1.7 x 10“^^ 
ergs“^ cm“^ arcsec“^) agree within the uncertainties^ 
The agreement between the Lya spectroscopic and nar¬ 
row band imaging surface brightnesses for both slit ori¬ 
entations confirm that our spectroscopic calibration pro¬ 
cedure is robust. 

We do not detect any extended emission in either the 
C IV or in the He II line, for either of the slit orientations. 

Note that one spatial dimension is set by the width of the slit, 
i.e. 1". 

We do not quote errors for these second set of measurements 
because there are significant systematics associated with the PSF 
subtraction in both imaging and spectroscopic data. 


































Dense Compact Gas Clumps in the CGM of UM287 


5 


To better visualize the presence of extended emission, we 
first subtract the PSF of the QSOs for each position an¬ 
gle (see middle rows in Figures and|^, and finally, we 
show in the upper rows the smoothed Xsmth maps. These 
smoothed maps are of g reat assistance in identifying fa int 
extended emission (see 


Hennawi & Prochaska 2013 for 


more details on the Psh' subtraction and the calcula- 
tion of the smoothed y-maps). The lack of compelling 
emission features in the PSF-subtracted smoothed maps 
confirm the absence of extended G IV and He II at our 
sensitivity limits in both slit orientations. 

As our goal is to measure line ratios between the Lya 
emission and the C IV and He II lines, we compute 
the surface brightness limits within the same aperture in 
which we calculated the Lya emission along the ’red’ slit, 
i.e. I"x20"and Av = 3000 km s“^. Because the com¬ 
panion quasar is much fainter than the UM287 quasar, 
the PSF subtraction along the ‘red’ slit does not suf¬ 
fer from systematics, whereas the large residuals in the 
left panel of Figure indicate that there are significant 
PSF subtraction systematics for the Lya emission in the 
‘blue’ slit covering the UM287 quasar. We have thus 
decided to focus on the line ratios obtained from the 
‘red’ slit, although the constraints we obtain from the 
‘blue’ slit are comparable. We find = 3.3 x 

10“^® ergs“^ cm“^ arcsec”^ and = 3.7x 10“^® 

ergs“^ cm“^ arcsec”^, respectively at the CIV and Hell 
locations. 

To better understand how well we can recover emission 
in the He II and C IV lines in comparison to the Lya, we 
visually estimate the detectability of extended emission 
in these lines by inserting fake sources as follows. First, 
we select the Lya emission above its local Icr limit along 
the ‘red’ slit, we smooth it and scale it to be 1, 2, 3, and 
5x at the location of the Hell and CIV line. Fi¬ 

nally, we add Poisson realizations of these scaled models 
into our 2-d PSF and sky-subtracted images. In Figure 0 
we show the x-maps for this test at the location of HelL 
This test suggests that we should be able to clearly de¬ 
tect extended emission on the same scale as the Lya line 
if the source is > 3x Thus, in the remainder of 

the paper we use 3cr (cr =SB^^^°) upper limits on the 
He II/Lya and C IV/Lya ratios. Given the values for the 
SBLya and the surface brightness limits at the location 
of tire CIV and Hell lines (within the I"x20"aperture 
and 3000 km s“^ velocity window for the red slit) we get 
(He II/Lya) 3 cr < 0.18 and (C IV/Lya) 3 cr < 0.16. Note 
that given the brighter Lya emission at the location of 
the ‘blue’ slit, the limits implied are about 2x lower than 
these quoted limits for the ‘red’ slit. 

It is important to note that we detect extended C IV 
emission around the faint companion quasar ‘b’ (see the 
smoothed maps in Figure]^. As this line is physically 
distinct from the UM287 nebula and essentially follows 
the extended Lya emission around the faint quasar (com¬ 
pare the smoothed maps for Lya and C IV), this suggests 
that we have detected the extended narrow emission line 
region (EELR) of this source. This kind of emission, pro¬ 
duced by the gas excited by an AGN on scales of tens 
of kpc, is u sually observed around low redshift z < 0.5 

2006 iHusemann et al. | |20I3 l 


the emission because, given the much smaller scales in 
play here, its accuracy depends on the PSF-subtraction. 
However note that this detection, near the limit of our 
sensitivity, clearly demonstrates that we could have de¬ 
tected faint extended emission in the C IV and He II 
lines within the Lya nebula itself if this emission were 
characterized by higher line ratios. 

Finally, we briefly comment on the nature of two other 
sources which fall within the ‘red’ slit, i.e. a Lyman 
Alpha emitter (LAE) (i.e. > 20A) and a contin¬ 

uum source (see Figure Indeed, this slit orientation 
was also chosen to confirm the presence of a LAE at 
about 350 kpc northward of ‘QSO b’, clearly visible in 
the na rrow-band image in Figure 2 of |Gantalupo et all 
(|20I4|) and in our Figure [11 Our LRIS data confirm the 
presence of a line emission from a LAE at a redshift 
z = 2.280 ± 0.002, which is consistent with the redshift 
of the UM287 quasar, within our uncertainties. We as¬ 
cribe this emission to the Lya line, and we compute a 
flux of F^ya = (9.2 ± 0.9) X 10“^® erg s“^ cm“^ (in 
an aperture of Av = 1400 km s“^, and 4 arcsec^), in 
agreement with F^ya = (8.4±0.4) x 10“^® erg s“^ cm“^ 
computed in an ape rture of 4 arcsec^ in the map of Can- 
talupo et al. (20141. We also serendipitously obtained a 
spectrum of a source at ~ 230 kpc from quasar ‘b’, which 
coincid es with a continuum so urces in our deep V-band 
image ( Cantalupo et al.||20I4 ). In our 2-d spectrum, we 
detect a faint continuuni associated with this source and 
an emission line at a wavelength of 5123A, which ap¬ 
pears at a velocity ^ 2750 km s“^ from the C IV line 
in the right panel of Figure However, given the low 
signal-to-noise ratio of the continuum, and the detection 
of a single emission line, we are unable to determine the 
redshift of this source. 

3.1. Kinematics of the Nebula 

With these slit spectra for two orientations, we can be¬ 
gin to study the kinematics of the Lya emission of the 
UM28'i^iant nebula. We first focus on the ‘red’ slit (see 
Figure^, which covers the companion quasar (‘QSO b’) 
and the extended Lya emission at a projected distance of 
100-160 kpc (~ 13"-19") from UM287 (‘QSO a’). We 
tested the kinematics of the detected emission by measur¬ 
ing the flux-weighted line centroid and the flux-weighted 
velocity dispersion (tr) around the centroid velocity in 2- 
pixels wide bins (^ 0.54") across the spatial slit direction. 
We then converted the velocity dispersion to a gaussian- 

~ 2.35cr. 


type-I (e.g. Stockton et ah 


and type-H quasars (e.g. Greene et al.|[2011), traced by 


[O III] and Balmer lines. We do not quote a value for 


equivalent FWHMgauss assuming FWHMg 
Note that, because of the resonant nature of the Lya 
emission, the line wi dth may be broadened by radiative 
transfer effects (e.g., Cantalupo et al. 2005) and repre¬ 
senting, thus, only an upper limit tor the thermal or kine- 
matical broadening. The extended emission has an aver¬ 
age FWHMgauss 500 km s“^ at a redshift of z = 2.279, 
which is centered on the systemic redshift of the UM287 
quasar. Although the emission appears coherent on this 
large scales, the gaussian FWHM calculated at each lo¬ 
cation ranges between ~ 370 km s“^ and ^ 600 km s“^, 
suggesting the need of higher resolution data to better 
characterize its width and shape. The line emission is 
red-shifted by ~ 750 km/s from quasar ‘b’. However note 
that our estimate for the redshift of quasar ‘b’ z = 2.275 
has a large 800 km s“^ error, because it is estimated 
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Figure 2. Two-dimensional spectra for the red slit shown in Figure^ plotted as x-uiaps following [Hennawi fe Prochask^ ( |2013[ |. In all 
panels, v = 0 km s“^ indicates the systemic redshift of the UM287 quasar, while the distance is computed from the companion quasar, 
i.e. ‘QSO b’. Bottom row: Xsky (sky-subtracted only) at the location of Lyo, He II, and C IV. Middle row: Xsky-i-PSF (sky and PSF 
subtracted) at the location of Lyo, He II, and C IV. Upper row: smoothed maps Xsmth after the PSF subtraction of the companion QSO 
(‘QSO b’ in Figure^. As expected, the extended Lyo emission is well visible in these panels up to 200 kpc from the companion QSO. 
Note also that within this slit we have a continuum source (source ‘C’ in Figure]^ at ~ 230 kpc, and a Lyman Alpha emitter (‘LAE’, also 
highlighted in Figureat ~ 350 kpc (see Section [j^for details). The blue box indicates the aperture used to compute the SBny ’ 
the limits on He Il/Lya and C IV/Lya line ratios, i.e. I"x20"and An = 3000 km s“^. 


and 


from broad rest-frame UV em ission lines which are p oor 
tracers of the systemic frame (Cantalupo et al.||2014). 

As for the ‘blue’ slit, statements about kinernatics are 
limited by the challenge of accurately subtracting the 
PSF of the bright UM287 quasar. Given our SB limit, we 
detect the Lya emission out to ~ 150 kpc. As expected 
from the narrow-band imaging, the Lya is stronger at 
this location in comparison with the other slit orienta¬ 
tion. In particular, the emission shows a peak at ^63 
kpc (~ 7.7") in agreement wi th the narrow-band dat a 
(see Figure or Figure 


FWHMgauss ~ 920 km s“^ and appears to vary more 
with distance along the slit. This larger width may arise 
from the fact we are probing smaller distances from the 
UM287 quasar than in the ‘red’ slit. 

Note that, at our spectral resolution (FWHM ^ 320 
km s“^), there is no evidence for “double-peaked” kine- 
matics characteristic of resonantly-trapped Lya (e.g. 
Cantalupo et al. 2005) along either slit. This may in- 


Cantalupo et al. 2014). 


At this second location, the Lya line appears broader 


dicate that resonant scattering of Lya photons does not 
play an important role in the Lya kinematics, however, 
data at a higher resolution are needed to confirm this 
conclusion. 
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Figure 3. Two-dimensional spectra for the blue slit shown in Figure^ plotted as x-uiaps following [Hennawi Prochask^ | |2013[ |. In 
all panels, v = 0 km s“^ indicates the systemic redshift of the UM287 quasar. The distance is also computed from the UM287 quasar, 
i.e. ‘QSO a’. Bottom row: Xsky (sky-subtracted only) at the location of Lyo, He II, and C IV. Middle row: Xsky-i-PSF (sky and PSF 
subtracted) at the location of Lyo, He II, and C IV. Upper row: smoothed maps Xsmth after the PSF subtraction of the UM287 QSO 
(‘QSO a’ in Figure [y. As expected, also along this slit we detect extended Lyo emission. Given our sensitivity limits, the Lyo line is 
detected up to ~15(jKpc from the UM287 QSO. Note that for such a brif^ QSO, it is difficult to cleanly subtract its PSF. The blue box 
indicates the aperture used to compute the SBnya as outlined in section O 


These estimates for the widths of Lya emission are 
comparable to the velocity widths observed in ab¬ 
sorption in the C GM surrounding z ~ . 2 quasars 


4. MODELING THE LYo, GIV AND HEH EMISSION 
AROUND UM287 


et al. 


(Av 500 k m s Prochaska & Hennawi 


2009 


Lau 


2015), pherapis suggesting that the kinematics 


As shown by jCantalupo et al.j (|2014|) 

IM; 


traced m emission are dominated by the motions of the 
gas as opposed to the effects of radiative transfer. Both 
the emission and absorption kinematics are comparable 
to the virial velocity ~ 300 km s“^ of the massi ve dark 
matter halo s hosting quasars (Mdm 10^^'® M^, [White 


et al. 2012), and thus appear consistent with gravita- 
tional motions. 


the extended 

Lycr emission nebula around UiVl287 can be explained 
by photoionization from the central quasar, imply¬ 
ing a large amount of cool (T ~ 10"^ K) gas, i.e. 
Me — 10^^ Mq. To further constrain the prop¬ 

erties of the gas in this huge nebula, in this section we 
exploit the si mple model for cool clouds in a quasar halo 
introduced by Hennawi & Prochaska (2013) and the con- 
se quent photoionization modeling~ procedure introduced 
by Arrigoni Battaia et al. (2014bI. Our main goal is to 
show how our line ratio constraints on G IV/Lya and 
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Figure 4. Illustration of detection significance of scaled models 
of the Lya emission at the location of the Hell line along the ‘red’ 
slit (Figure]^. The synthetic sources corresponds to 1, 2, 3, and 
5x bottom panel shows the Xsky (sky-subtracted 

only) map, while the upper panel shows the smoothed map. This 
figure suggests that we should be able to clearly detect extended 
emission > 3cr on the scale of the Lya line. 

He II/Lya can be used to constrain the physical proper¬ 
ties of the gas in the UM287 nebula, such as the volume 
density (u-h), column density (Xh), and gas metallicity 
(Z). 


We reiter ate that as in our previous work (Cantalupo^ 
et al.|[2M4 ), model the Lya emission alone cannot break 
the degeneracy between the dumpiness or density of the 
gas, and the total gas mass. In the next sections we show 
how information on additional lines (in particular He II) 
can constrain the density of the emitting gas and thus 
break this degeneracy. 

4.1. Photoionization Modeling 

In the following, we briefly ou tline the simple model 
for co ol halo gas introduced by Hennawi & Prochaska 
(2013) for the case of UM287. We assume a simple pic- 


ture where UM287 has a spherical halo populated with 
spherical clouds of cool gas (T ^ 10^ K) at a single uni¬ 
form hydrogen volume density nn, and uniformly dis¬ 
tributed throughout the halo. We model a scale length 
of i? = 160 kpc from the central quasar, which approx¬ 
imately corresponds to the distance probed by the ‘red’ 


slit, and represents the expected virial radius for a dark 
matter halo hosting a quasar at this redshift. In this 
configuration, the spatial distribution of the gas is com¬ 
pletely specified by nn, 7?, the hydrogen column density 
Xh, and the cloud covering factor fc- 
Note that the tot al mass of cool gas in our simp le model 
can be written as (Hennawi fc Prochaska|[^013[): 


M, = TTR^fcNu 


= 2.7 X 10 


10 


A 

/ R 


( 2 ) 


Nh 


\ 160 kpc/ \ 10^®-^ cm“2 / \^ 1.0 


fc 


Mp 


where mp is the mass of the proton and X is the hydrogen 
mass fraction. 

In this simple model, the Lya SB is determined by 
simple relations which depend only on nn, Xh, fc, and 
the luminosity of the QSO at the Lyman limit (L;yLL)(see 
Hennawi &: Prochaska 2013 for details). To build intu- 
ition, it is useful to consider two limiting regimes for 
the recombination emission, for which the clouds are op¬ 
tically thin (Xhi ^ 10 ^^'^ cm“^) and optically thick 
(Ahi ^ 10^^'^ cm“^) to the Lyman continuum photons, 
where Xhi is the neutral column density of a single spher¬ 
ical cloud. We argue below, that given the luminosity of 
the UM287 quasar, the optically thick case is however 
unrealistic. 

- Optically thin to the ionizing radiation: 

SBI-a, (l + T) ,3, 

where Tythin = 0.42 is the fraction of recombinations 
which result in a Lya photon in t he optically thin 
limit (Osterbrock & Ferland|2006 ), h is the Planck 
constant, t'Lya is the frequency of the Lya line, 
aA = 4.18 X 10“^^ cm“^ s“^ is the c ase A recombi- 
nation coefficient at T = 10, 000K ( Osterbrock fc| 


Ferland 


2006 


and X = 0.76 and Y = 0.24 are 


implied by Big Bang Nucleosynthesis (jBoesgaard 

& Steigman 119851 

Izotov et al.| 

1999 

locco et al. 

2009 Planck Colla 

roration et a 

I 20 T 4 

). 


Optically thick to the ionizing radiation 
■/C$LL [R/Vs] 


OTiithick _ ^thick^^Lya 
“ 47r(l + z)4' 


( 4 ) 


where r?thick = 0.66 is the fraction of ionizing pho¬ 
tons converted into Lya photons, and where 4 >ll 
([ phot s“^ cm“^]) is the ionizing photon number 
flux. 

Thus, in the optically thick case the Lya surface bright¬ 
ness scales with the luminosity of the central source, 

Note that this equation hides a dependence on temperature 
through the recombination coefficient qa* which usually is ne¬ 
glected, but that can be important, i.e qa is decreasing by a factor 
of ^6 froni T = 10^ K to T = 10^ K (CHIANTI database, 

|et al.|1997[[Landi et al.|2013[ |. 
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oc fcLy^^, while in the optically thin regime the 
SB does not depend on oc /c’T'H-^h, pro¬ 

vided the AGN is bright enough to keep the gas in the 
halo ionized enough to be optically thin. 

We now argue that the Lya emitting gas is unlikely to 
be optically thick Ahi > cm“^. Equationsand 

can be combined to express the SB in terms of , the 
luminosity at the Lyman edge. To compute this lumi¬ 
nosity, we assume that the quasar spectral energy distri¬ 
bution obeys the power-law form L^, = , 

blueward of t'LL and adopt a slo pe of ottv = —TT consis¬ 
tent with the measurements of Lusso et al.| (|2015 ). The 
quasar ionizing luminosity is then parameterized by Tj/ll ? 
the specific luminosity at the Lyman edg^^ We deter¬ 


mine the no rmalization Lj/ll by integrating the Lusso 


et al. (2015) composite spectrum against the SDSS hlter 


curve, and choosing the amplitude to give the correct i- 
band magnitude of the UM287 quasar (i-mag= 17.28), 
which gives a value of = 5.4 x 10^^ erg s“^ Hz“^. 

Substituting this value of for UM287 into equa¬ 
tion]^ we thus obtain 


SB^)';))‘' = 8.8 X 10-1® 


1 


3.279 


-4 


-2 


r5)(T6,^)' 


X 


Ly 


This value is over two order of magnitude larger than the 
observed SB value of the Lya emission at 160 kpc from 
UM287. Even if we consider a larger radius, R = 250 kpc, 
in order to get the observed SBi^ya we would need a very 
low covering factor, i.e. fc ^ 0.02. Such a small covering 
factor would be strictly at odds with the observed smooth 
morphology of the diffuse nebula as seen in Figured We 
directly test this assumption as follows. We ran^mly 
populate an area comparable to the extent of the Lya 
nebula with point sources such that fc = 0.1 — 1.0, and 
we convolve the images with a Gaussian kernel with a 
FWHM equal to our seeing value, in order to mimic the 
effect of seeing in the observations. We find that the 
smooth morphology observed cannot be reproduced by 
images with fc < 0.5, as they appear too clumpy. Thus, 
the smooth morphology of the emission in the Lya nebula 
implies a covering factor of fc ^ 0.5. 

In the following sections we construct photoionization 
models for a grid of parameters governing the physical 
properties of the gas to estimate the expected He II and 
C IV emission. Following the discussion here, we shall 
see that the models which reproduce the observed Lya 
SB will be optically thin, because given the high covering 
factor optically thick models would be too bright. 


4.2. The Impact of Resonant Scattering 

It is important to stress at this point that the Lya 
photons should be subject to substantial resonant scat¬ 
tering under most of the astrophysical conditi ons, given 


the large optica l depth at line center (see e.g. Gould & 


Weinberg 1996). Thus, typically, a Lya photon experi- 
ences a large numbers of scattering before escaping the 
system in which it is produced. This process thus leads 


We describe in detail t he as sumed quasar spectral-energy dis¬ 
tribution (SED) in Section j]4.3| 


to double-peaked emission line profiles as Lya photons 


be able to escape the system (e.g. 

Neufeld 1990 Gould 

& Weinberg||1996l ICantalupo et al. 

|2UU5t IDijkstra et al. 

2006b| 

Verhamme et al.||2006). Altl 

lOugh our models are 
i.e. to ionizing pho- 

optical 

ly thin at the Lyman limit. 


tons, for the model parameters required to reproduce the 
SB of the emission, they will almost always be optically 
thick to the Lya transition (i.e. iVni > 10^"^ cm”^). 
Hence one should be concerned about the resonant scat¬ 
tering of Lya photons produced by the central quasar 
itself. However, radiative transfer simulations of radi¬ 
ation from the UM287 quasar through a simulated gas 
distribution have shown that the scattered Lya line pho¬ 
tons from the quasar do not contribute significantly to 
the Lya surface brightness of the nebu la on large scales. 


i.e. > 100 kpc (Gantalupo et al. 2014). This is because 
the resonant scattering process results in very efficient 
diffusion in velocity space, such that the vast majority of 
resonantly scattered photons produced by the quasar it¬ 
self escape the system at very small scales < 10 kp c, and 


stra et al.|2006b[[VeHiamme et al.j2006[[C 

2005|). For this reason, based on the result 


hence do not propagate at larger distances (e.g. Dijk- 

antalupo et al. 

results of the radia¬ 


tive transfer simulations of Gantalupo et al. (2014), we 
do not model the contribution of resonant scattering of 
the quasar photons to the Lya emission. Similar consid¬ 
erations also apply to the resonant C IV line, however we 
note that resonant scattering of C IV is expected to be 
much less efficient, because the much lower abundance of 
metals imply the gas in the nebula is much less likely to 
be optically thick to C IV. 

To avoid a contribution to the Lya and G IV emission 
from scattering of photons from the QSO we have thus 
masked both lines in our assumed input quasar spec¬ 
trum. Note that with this approach we do not neglect 
the ‘scattered’ Lya photons arising from the diffuse con¬ 
tinuum produced by the gas itself, which however turn 

out to be insignificant E3 


4.3. Modeling the UM287 Quasar SED 

We assume that the spectral energy distribution (SED) 
of UM287 has the form shown in Figure]^ As we do not 
have complete coverage of the spectrum of this quasar, 
we adopt the following assumptions to model the full 
SED. Given the ionization energies for the species of in¬ 
terest to us in this work, i.e. 1 Ryd=13.6 eV for Hydro¬ 
gen, 4 Ryd=54.4 eV for He II, and 64.5 eV for C IV, we 
have decided to stick to power-law approximations above 



us to relate the i-band magnitude of the UM287 quasar 
to the specific luminosity at the Lyman limit ■ For 


Note that this value depends on the broadening of the line 
due to turbulent motions of the clouds. Given current estimates 
of typical equivalen t widths of optically thi ck absorbers in quasar 
spectra, i.e. ~ lA ( jProchaska et ah |2013b| | , in our calculation we 
consider turbulent motions of 30 km s~^. However, note that our 
results are not sensitive to this parameter. 
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energies greater than one Rydberg, we assume a power 
law form Ly = and adopt a slope of 


qttv = —1.7. consistent with the measurements of Lusso 


et al. (20151, while in the Appendix we test also the 


cases for cteuv = and —2.3. We determine the 

normalization by integrating the Lusso et al. (20151 
composite spectrum against the SDSS hlter curve, and 
choosing the amplitude to give the correct i-band mag¬ 
nitude of the UM287 quasar (i.e. *=17.28), which gives 
a value of = 5.4 x 10^^ erg s“^ Hz“^. We extend 
this UV power law to an energy of 30 Rydberg, at which 
point a slightly different power law is chosen a = —1.65, 
such that we obtain the correct value for the specific lu¬ 
minosity at 2 keV L^(2keV) implied by measurements 
of aox, defined to be Ly(2keV)/Li/(2500 A) = (j^ 2 kev/ 
^adopt the value aox = -1.5 measured 

by Strateva et al. (2005) for SDSS quasars. An X-ray 
slope of ax = —1, which is flat in vfi, is adopted in 
the interval of 2-100 keV, and above 100 keV, we adopt 
a hard X-ray slope of cthx = —2. For the rest-frame op¬ 
tical to mid-lR part of the SEP, we spli c e together the 


composite s pectr a of Lusso et al .| ( 201 5), Vanden Berk 

2006|^. 'I'hese assump- 


et al. (2001), and Richards et al. 
tions about the SEP are essentially tbe standa rd ones 


used in ph otoionization modeling of AGN (e.g. Baskin 
et al.|2014 ). Summarizing, given the lack of information. 


tor energies greater than one Rydberg we parametrized 
the SED of the UM287 quasar with a series of power-laws 
as 


/i/ oc 


jyOCBW 

hv >\ Ryd 


if 30 Ryd <hv <2 keV 


if 2 keV <hv< 100 keV 


if hv > 100 keV. 


(7) 


line, we also consider variations in Z. Thus, we run a uni¬ 
form grid of models with this wide range of parameters: 

- riH = 10“^ to 10^ cm“^ (steps of 0.2 dex); 

- Ah = 10 ^® to 10^2 (steps of 0.2 dex); 

- Z = 10“® Zq to Zq (steps 0.2 dex). 


Note that by exploring this large parameter range, some 
of the models that we consider result in clouds opti¬ 
cally thick at the Ly man limit, but as explained in 
the previous Section §4.1[ these parameter combina¬ 
tions result in nebulae which are too bright and thus 
inconsistent with the observed Lya surface brightness. 
In what follows, we only consider the models which 
closely reproduce the observed Lya surface brightness, 
i.e. 5.5 X 10“^® ergs“^ cm“2 arcsec“2 < SBLya < 8.5 x 

10 “^® ergs“^ cm“2 arcsec“2. 

Photoionization models are self-similar in the ioniza¬ 
tion parameter U = which is the ratio of the num¬ 
ber density of ionizing photons to hydrogen atoms. As 
the luminosity of the central QSO is known, the variation 
in the ionization parameter U results from the variation 
of the volume number density rin for the models in our 
grid. The range of ionization parameters that we cover 
is comparable to those in previous analysis of photoion¬ 
ization around AG Ns, e.g. in the case of the narrow line 
regions (NLR; e.g. Groves et al.|2004) and in the case of 


extended e mission line regions (EELR; e.g. Humphrey 
et al.|[2008|. Finally, we emphasize that once we hx the 


source luminosity and define the ionizing spectrum, the 
line ratios we consider are described by two model param¬ 
eters, namely the density nn of the gas and its metallicity 
Z. We will see this explicitly in the next section. 


4.4. Input Parameters to Cloudy 

Having established our assumptions on the UM287 
SED, and on the resonant scattering, we now explain 
how we choose the range of our model parameter grid. 
We perform our calculations with the Gl qudy photoion¬ 
ization code (vlO.Ol), last described by Ferland et al. 
(2013). Because the emitting clouds are expected to be 
much smaller than their distance r ^ i?vir = 160 kpc 
from the central ionizing source, we assume a standard 
plane-parallel geometry for the emitting clouds illumi¬ 
nated by the distant quasar. In order to keep the models 
as simple as possible, and because we are primarily in¬ 
terested in understanding how photoionization together 
with the observed line ratios can constrain the physi¬ 
cal properties of the gas (i.e. nn and Ah), without 
resorting to extreme parameter combinations, we pro¬ 
ceed as follows. We focus on reproducing the SBLya ^ 
7 X 10“^® ergs“^ cm“2 arcsec”^ at 160 kpc distance from 
the UM287 quasar which is basically the distance probed 
by the ‘red’ slil[^ In particular, eqn. ^ implies that 
a certain combination of Ah and nn are thus required. 
Further, given the dependence on metallicity {Z) of the 
G IV and He II lines, and of the gas temperature which 
determine the amount of collisional excitation in the Lya 


Note that we have decided to model a single distance from 
the UM287 quasar. The sensitivity of our results to this simple 
assumption is discussed in Section [tI 


5. MODELS VS OBSERVATIONS 

As we discuss in Section our LRIS observations pro¬ 
vide upper limits on the C IV/Lya and He II/Lya ratios, 
i.e. (G IV/Lya) 3 CT < 0.16 and (He II/Lya) 3 cr < 0.18. On 
the other hand, each photoionization model in our grid 
predicts these line ratios, and Figure shows the trajec¬ 
tory of these models in the He II/Lya vs C IV/Lya plane. 
The region allowed given our observational constraints on 
the line ratios is indicated by the green shaded area. We 
remind the reader that we select only the models which 
produce the observed Lya emission of SBpya ~ 7 x 10“^® 
ergs“^ cm“2 arcsec“2^ which to lowest order requires a 
combination of Ah and nn as shown by eqn. <§■ Since 
the luminosity of the central source is known, these mod¬ 
els can be thought to be parametrized by either nn or 
the ionization parameter U, as shown by the color coding 
on the color-bar. In the same plot we show trajectories 
for different metallicities Z = 1, 0.1, 0.01, 10~^ Zq. 

W e now reconsider the covering factor. We argued in 
§4.1| that based on the morphology of the nebula, the 
covering factor need to be fc ^ 0.5, and that optically 
thick gas clouds would tend to overproduce the Lya SB 
for such high covering factors. Our models provide a 
conhrmation of this behavior. For a covering factor of 
fc = 1.0 a large number of models are available, whereas 
if we lower the covering factor to fc = 0.3, we find that 
only two models in our extensive model grid can satisfy 
the Lya SB of the nebula. This results because as we 
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Figure 5. Spectral energy distribution (SED) of UM287 used as incident radi ation field in our m odeling. The black solid line indicate our 
fiducial input spectrum characterized by a slope in the EUV of ckeuv = —1.7 l |Lusso et al. |2015| l. The right panel is a zoomed version of 
the box highlighted in the left panel. Note the agreement between the composite spectrum used as input in our calculation and the SDSS 
spectrum of UM287 (green solid line). To prevent a contribution from resonantly scattered photons, we mask the emission from the line 
center of Lya and C IV. 


decrease fc, assuming the gas is optically thin, eqn. ^ 
indicates we must correspondingly increase the product 
Vh^ih by l//c in order to match the observed Lya SB. 
However, note that the neutral fraction also scales with 
this product xm oc Vh^ih such that for low enough val¬ 
ues of fc increasing iVnriH would result in self-shiel ding 
clouds that are optically thick. We already argued in §4.1| 
that if the clouds are optically thick the covering factor 
must be much lower fc — 0.02, which is ruled out by the 
diffuse morphology of the nebula. Hence our constraint 
on the covering factor fc ^ 0.5 can also be motivated by 
the simple fact that gas distributions with lower cover¬ 
ing factors would over-produce the Lya SB. Henceforth, 
for simplicity, we assume a covering factor of fc = 1-0 
throughout this work, but in ^we test the sensitivity of 
our results to this assumption. 

The gray symbols in Figure also show a compila¬ 
tion of measurements of the He II/Lya and C IV/Lya 
line ratios from the literature for other giant Lya neb- 


(2014b 

). Specifically, we show measurements or u 

pper 

limits 

for the two line ratios for seven Lya blobs 1 

De;y 

et al. 

20051 Prescott et al. 2009 |2013[ |Arrigoni 

Bat- 


5d high redshitt rac 


Villar-Martm et al. 


io ga laxies (Humphrey et al. 2006 


2007). Note that we show measure¬ 


ments from the literature in Figure!^ for reference, but 
these measurements cannot be directly compared to our 
observations or our models for several reasons. First, the 
emission arising from the narrow line region of the cen¬ 
tral obscured AGN is typically included for the HzRGs, 


From the sample of|Arrigoni Battaia et al.| |2014bJ . We decide 
to plot only the upper limits ot UABi and LAB2, which set the 
tighter constraints for that sample. 


contaminating the line ratios for the nebulae. In addi¬ 
tion, the central source UV luminosities are unknown for 
both LABs and HzRGs, and thus they cannot be directly 
compared to our models, which assume a centr al source 
luminosity. See |Arrigoni Battaia et al.| 2014b| and ref¬ 
erences therein for a discussion on this dataset and its 
caveats. 

The trajectory of our optically thin models through the 
He II/Lya and G IV/Lya diagram can be understood as 
follows. We first focus on the curve for Z = Zq and fol¬ 
low it from low to high U (i.e. from high to low volume 
density tt-h). First consider the trend of the He II/Lya 
ratio. He II is a recombination line and thus, once the 
density is fixed, its emission depends basically on the 
fraction of Helium that is doubly ionized. For this rea¬ 
son, the He II/Lya ratio is increasing from logG ^ —3 
and ‘saturates’, reaching a peak at a value of ^ 0.34 
which is set by atomic physics and in particular by the 
ratio of the recombination coefficients of Lya and He II. 
Indeed, if we neglect the contribution of collisional ex¬ 
citation to the Lya line emission, which is a reasonable 
assumption near solar metallicity, then both the He II 
and Lya are produced primarily by recombination and 
the recombination emissivity can be written as 


/line — /v 


elem 


^^line 

47r 


^e^ion^line (^) ? 


( 8 ) 


where riion is the volume density of Heand H+ for 
the case of Hell and Lya, respectively. Here afj®g(T) is 
the temperature dependent recombination coefficient for 
Hell or Lya, and the factor = 3/cIVeiem/(4i?neiem) 
takes into account that t he emitting clouds fill only a 
fraction of the volume (see Hennawi fc Prochaska|2013). 


Thus, once the Helium is completely doubly ionized, i.e. 
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Figure 6. Hell/Lya versus CIV/Lya log-log plot. Our upper 
limits on the Hell/Lyo and CIV/Lya ratios are compared with 
the Cloudy photoionization models that reproduce the observed 
SBLya ~ 7 X 10“^® ergs“^ cm~^ arcsec“^. For clarity, we plot 
only the models with Z = 0.001,0.01,0.1, and IZq. The models 
are color coded following the ionization parameter U, and thus 
the volume density nn (see color bar on the right). The green 
shaded area represents the region defined by the upper limits of 
the UM287 nebula. Note that these upper limits favor models 
with riH ^ 3 cm“^, A^h £10^*^ cm“^, and logU < —1.5. This is 
even more clear in Figure}^ 


np ~ riH and nHe++ ~ (V/2V)nH, the ratio between the 
two lines is given by the relation 


^= 0.34 

jhya 


<„(20,000K) \ 

1.15 X 10“i2 cm^ s“^ / 


/ al%i20,000K) \ 

( 2.51 X 10“i3 cm^ y ’ 


( 9 ) 


Note that eqn. ® depends slightly on temperature, with 
a decrease of the ratio at higher temperatures. Be¬ 
fore reaching this maximum line ratio, He II/Lya is 
lower because Helium is not completely ionized, and is 
roughly given by Hell/Lya ~ a;He++ x (jHeii/jLya)max, 
where a;He++ is the fraction of doubly ionized Helium. 
As stated above, this simple argument does not take into 
account collisional excitation of Lya. In particular, at 
lower metallicities when metal line coolants are lacking, 
the temperature of the nebula is increased, and collision- 
ally excited Lya, which is extremely sensitive to tem¬ 
perature, becomes an important coolant, boosting the 
Lya emission over the pure recombination value. Thus 
metallicity variations result in a change of the level of 
the asymptotic Hell/Lya ratio as seen in Figure]^ 

Our photoionization models indicate that the C IV 
emission line is an important coolant and is powered pri¬ 
marily by collisional excitation. The efficiency of C IV 
as a coolant depends on the amount of Carbon in the 
ionic state. For this reason, the C IV/Lya ratio is 
increasing from logH ~ —3, reaches a peak due to a max¬ 
imum in the fraction, and lowers again at higher U 
where Carbon is excited to yet higher ionization states, 
e.g. C V. For example, for the Z = 0.1 Zq models, 


the C IV/Lya ratio peaks at log U = —1.4 and then 
decreases at higher U. Given that C IV is a coolant, 
the strength of its emission depends on the metallicity of 
the gas. Indeed, for metallicities lower than solar, C IV 
becomes a sub-dominant coolant with respect to colli- 
sionally excited Lya (and for very low metallicity, e.g. 
Z = 10~^Zq, also to He Lya), and its emission becomes 
metallicity dependent as can be seen in Figure 

At lower metallicities the Lya line becomes an impor¬ 
tant coolant. For the Z = O.OOIZq grid, the collisional 
contribution to Lya has an average value of ~ 40%, 
while it decreases to ^ 37%, ^ 25%, ~ 1% for the 
Z = 0.01,0.1,1^0 cases, respectively. Given that the 
strength of the collisionally excited Lya emission in¬ 
creases with density along each model trajectory, this 
slightly dilutes the aforementioned trends in the He II 
and C IV line emission. Specifically, the density depen¬ 
dence of collisionally excited Lya emission moves the line 
ratios to lower values for logt/ > —1.5, which would 
otherwise asymptote at the expected He II/Lya ratio 
in eqn. 0- Thus the effect of collisionally excited Lya 
emission tend to mask the ‘saturation’ of the He II/Lya 
ratio due to recombination effects alone, and results in a 
continuous increase of He II/Lya with U. 

Overall, Figure [^illustrates that our simple photoion¬ 
ization models can accommodate the constraints implied 
by our observed upper limits on the He II/Lya and 
C IV/Lya ratios of UM287. In particular, our non¬ 
detections are satisfied (green shaded region) for models 
with high volume densities nn and low metallicities Z. 
These constraints can be more easily visualized in Fig¬ 
ure where we show the allowed regions in the nj^-Z 
plane implied by our limits on the He II/Lya (panel ‘a’) 
and C IV/Lya ratios (panel ‘b’). Specifically, in these 
panels the solid black line indicate the upper limits in 
the case of the UM287 nebula, i.e. He II/Lya < 0.18 
(or log(He II/Lya)< —0.74), and C IV/Lya = 0.16 (or 
log(C IV/Lya)< —0.79), while the arrows indicate the 
region of the parameter space that is allowed. It is ev¬ 
ident that our limits on the extended emission in the 
He II/Lya ratio give us stronger constraints than those 
from the C IV/Lya ratio. The He II/Lya ratio provides 
a constraint on the volume density which is metallicity 
dependent, however even if we assume a logj^Q Z ~ —2—3, 
which are the lowest possible values com parable to the 


backg round metallicity of the IGM (e.g. Schaye et al. 


20031, we obtain a conservative lower limit on the vol- 
ume density of rin ^ 3 cm“^. 

Given this constraint on tih, and the fact that we know 
the Lya emission level, which in turns approximately 
scales as uhA^h (see eqn. ®), we can use our lower limit 
on riH to place an upper limit on Ah or equivalently on 
the total cool gas mass because it scales as /c Ah once 
the radius is fixed (see eqn. (^). Panel ‘c’ of Figure 
shows that our limit on the He II/Lya ratio combinea 
with the total SBLya implies the emitting clouds have 
column densities Ah < 10^° cm“^. Thus, if we assume 
that the physical properties of the slab modeled at 160 
kpc are representative of the whole nebula, we can com¬ 
pute a rough estimate for the total cool gas mass. With 
this strong assumption, that nn ^ 3 cm“^ is valid over 
the entire area of the nebula , i.e. 911 arcsec^ (from the 
2o- isophote of the Lya map; Gantalupo et al.|[2014|, we 
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Figure 7. Constraints on the physical parameters of the gas clouds from our photoionization models that reproduce the observed 
SBLya ~ 7 X 10“^^ ergs~^ cm“^ arcsec“^ in the case of an input spectrum with oeuv = —1-7. Given the known luminosity of the central 
source, the assumed SED, and the fixed SBLyo:? our models can be thought to be parametrized by only two parameters, namely nu and 
Z. Panel ‘a’: map of the He Il/Lya ratio in the nn-Z plane. The black solid line indicate our 3cr upper limit He Il/Lyo < 0.18 (i.e. 
log(He H/Lya)< —0.74). Panel ‘b’: map of the C IV/Lyo; ratio in the nn-Z plane. The black solid line indicate our 3 (j upper limit 
C IV/Lya < 0.16 (i.e. log(He Il/LyQ:)< —0.79). Note that the constraints from the He Il/Lya ratio are stronger. Panel ‘c’: map of the 
He Il/Lya ratio in the Nn-Z plane. The black solid line indicate our 3cr upper limit. Models with A^h $ 10^^ cm“^ are selected. Panel 
‘d’: map of the He Il/Lya ratio in the R-Z plane. The black solid line indicate our 3a upper limit. Note that really small cloud sizes are 
favored, i.e. R <20 pc. 
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then deduce that -/Vh < 10 ^° cm“^ over this same area, 
and hence the total cool gas mass is Me ;< 6.4 x 10^° M©. 

Further, by combining the lower limit on volume den¬ 
sity riH and upper limit on column density A^h, we can 
also obtain an upper limit on the sizes of the emitting 
clouds defined as R = iVn/nn- Panel ‘d’ in Figure Jb 
shows that this upper limit is constrained to be i? < ^ 
pc. Assuming a unit covering factor fc = 1.0, this con¬ 
straint on cloud sizes implies > 53,500 clouds per square 
arcsec on the sky, and each cloud should have a cool 
gas mass < 1.3 x 10^ Mq. Assuming these clouds 
have the same properties throughout the whole nebula, 
we find that > 4.9 x 10^ clouds are needed to cover the 
extent of the Lya emission (~ 911 arcsec^ 

The foregoing discussion indicates that we are able 
to break the degeneracy between the volume density of 
the gas rtH and the to tal cool gas mass presented in 


Cantalupo et al. (2014). As a reminder, this degener- 
acy arises because the Lya surface brightness scales as 
SBLya oc nniVn, whereas the total cool gas mass is given 
by Me oc A^h- Thus observations of the Lya alon e can¬ 
not independ e ntly d etermine the cool gas mass. |Can-| 


talupo et al. (2014) modeled the Lya emission in the 
UA1287 nebula m a way that differs from our simple 
model of cool clouds in the quasar CGM. Specifically, 
they used the gas distribution in a massive dark matter 
halo M = 10^^ ® Mq meant to represent the quasar host, 
and carried out ionizing and Lya radiative transfer sim¬ 
ulations under the assumption the gas is highly ionized 
by a quasar with the same luminosity as UM287, and 
the extended Lya emission is dominated by recombina¬ 
tions, similarly to our simpler Cloudy model^^ Under 
these assumptions, they are not able to reproduce the ob¬ 
served Lya surface brightness of the nebula. This arises 
because only ~ 15% of the total gas in the simulated 
halo is cool enough to emit Lya recombination radiation 
(r < 5 X 10^ K), because the vast majority of the baryons 
in the halo have been shock-heated to the virial temper¬ 
ature of the halo, i.e. T ~ 10^ K. Even if they assume all 
of the gas in the simulated halo can produce t he Lyg line 
(Myas ~ 10 ^^'^ Mq for the dark matter halo; Cantalupo 


et al. 2014), the surface brightn ess of the resulting neb¬ 


ula is still too faint. As a result, Cantalupo et al. (2014) 
postulated that the emission in the simulated halo must 
be boosted by a clumping factor C = (ny)/(nH)^, which 
represents the impact of clumps of cool gas which are 
not resolved by the simulation. They then determined 
the scaling relation between the simulated Lya emission 
and the column density of the si mulated gas distribution . 


i.e. 


A^h oc SBLyo,^/^ C 1 ^ ( [Cantalupo et al. 


as expected for recombination radiation. iNote that ac- 


2014), 


We quote a lower limit on the number of clouds per arcsec^ 
because we calculate this value without taking into account the 
possible overlap of clouds along the line of sight, and also because 
we use the maximum radius allowed by our constraints. In other 
words, we simply estimate the number of clouds with radius R = 20 
pc needed to cover the area of a square arcsec on the sky at the 
systemic redshift of the UM287 quasar. 

Although note that our Cloudy models treat collisionally ex¬ 
ited Lya emission properly, whereas this effect cann ot be properly 
modele d via the method i n Can talupo et al.l 1 2014|l. 

In [Cantalupo et al.| |2 ll4| this relation is quoted as A^hh oc 
Rt 1/^^ but Wh ~ AIhii in this simulated case where the 


SBlv 


gas is highly ionized. 


cordingly, SBLya^^^ oc Ah oc Me and one sees 
that this is identical to the scaling implied by eqn. (j^, 
SBLyct oc NiiUii, if one identihes nn with Our 

simple cloud model adopts a single density for all the 
clouds tih, whereas in the clumping picture, there could 
be a range of densities present, but the emission is dom- 
inated by gas w i th np — (^h). In this context, 

Cantalupo et al. (2014) inferred that if C = 1, the high 
observed SBLya, implies very high column densities up 
to Ah « 10 ^^ cm“^ corresponding to cool gas masses 
Me = 10^^ Mq, in excess of the baryon budget of the 
simulation. More generally, in the presence of clumping 
this constraint becomes Me = 10^^ Mq. 

By introducing the constraint on the volume density 
3 cm“^ using the He II/Lya ratio, our work (i) 
breaks the degeneracy between density nn (or equiva¬ 
lently C) and total column density Ah (or equivalently 
Me), (ii) allows us to then constrain the total cool gas 
mass Me < 6.4 x 10^° Mq without making any assump¬ 
tions about the quasar host halo mass, and (iii) de¬ 
mands the existence of a population of extremely com¬ 
pact (i? < 20 pc) dense clouds in the CGM/IGM. The 
ISM-like densities and extremely small sizes of these 
clouds clearly indicate that they would be unresolved by 
current cosmological hydrodynam ical simulations, given 
their resolution on galactic scales (iFumagalli et al. 112014 

-- lion elll.yOi'.l Iii 


Faucher-Giguere et al.|2014{|CrighTOn et al.|2015[|lNeisoii| 


et al. 2015). Indeed, our measurements wo uld imply a 


ciumpingme tor C > 200 for the simulation of Cantalupo 
et al.| (2014), in agreement with the value they required 


in order to reproduce the observed Lya from their simu¬ 
lated halo. 


5.1. Constraints from Absorption Lines 

A source lying in the background of the UM287 neb¬ 
ula that pierces the gas at an impact parameter of ~ 160 
kpc may also exhibit absorption from high-ion UV tran¬ 
sitions like C IV and N V, which can be constrained 
from absorption spectroscopy. In Figure we show a 
map for the column density of the C IV and N V ionic 
states {New, -^nv) for our model grid that reproduces 
the observed SBlyq ~ 7 x 10“^® ergs“^ cm“^ arcsec”^. 
Given our non-detection of He II emission, our upper 
limits on the He II/Lya ratios (indicated by the black 
solid line in both panels), imply Aqiv 10^® ® cm“^ and 
Anv ^ 10^® ° cm“^, respectively. The quasar UM287 re¬ 
sides at the center of the nebula, and our narrow band 
image indicates it is surrounded by Lya emitting gas. It 
is thus natural to assume that the UM287 quasar pierces 
the nebular gas over a range of radial distance^^ Thus 
a non-detection of absorption in these transitions places 
further constraints on the physical state of the absorbing 
gas in the nebula. 

To this end, we examined the high signal-to-noise 
S/N ~ 70 pix“® SDSS spectrum of the UM287 quasar, 
which has a resolution of i? ~ 2000. We hnd no evidence 
for any metal-line absorption within a ~ 2000 km s“® 
window of the quasar systemic redshift coincident with 
the velocity of the Lya emitting nebula (see Figure 

This would not be the case if the emitting gas is all behind the 
quasar. Given that the quasar shines towards us and contemporary 
on the gas, this configuration seems unlikely. 
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Figure 8. Top Panel: Map of the C IV column density A^civ 
the riYi-Z plane built from our photoionization models that repro¬ 
duce the observed SBLyc^ ~ 7x 10“^^ ergs“^ cm“^ arcsec“^ in the 
case of an input spectrum with oeuv = —1-7. The black solid line 
indicate our 3 (j upper limit in the He Il/Lya ratio, while the gray 
dashed line indicate our limit of A^civ < cm~^ implied by 

the absence of absorption at the resolution of the SDSS spectrum of 
UM287. Bottom Panel: Map of the N V column density A^nv 
the nu-Z plane built from our photoionization models that repro¬ 
duce the observed SBLya ~ 7x 10“^® ergs“^ cm~^ arcsec”^ in the 
case of an input spectrum with oeuv = —1.7. The black solid line 
indicate our 3 (t upper limit in the He Il/Lya ratio, while the gray 
dashed line indicate our limit of A^nv < cm“^ implied by 

the absence of absorption at the resolution of the SDSS spectrum 
of UM287. The spectroscopic constraints for both species imply 
that the gas along the sightline, if present, is in a similar state 
as the observed nebula, being illuminated by the bright quasar as 
well. 


[^, implying IVciv < cm“^ (EWciv < 15 mA), 

and -/Vnv < cm“^. These limits constrain the 

amount of gas in these ionic states intercepted by the 
quasar at all distances, but in particular at ~ 160 kpc, 
where we conducted our detailed modeling of the emis¬ 
sion. As such, directly analogous to our constraints from 
the emission line ratios, we can similarly determine the 
constraints in the n^i-Z plane from the non-detections of 
C IV and N V absorption, which are shown as the gray 
dashed lines in Figurej^ As expected these metal absorp¬ 
tion constraints depend sensitively on the enrichment of 
the gas, but the region of the nn-Z plane required by 
our non-detections are consistent with that required by 
our He Il/Lya emission constraint. Specifically, for log 
Z > —2.3, the absence of absorption provides a compa¬ 
rable lower-limit on the density as the non-detection of 
emission, whereas at lower metallicities the absorption 
constraint allows lower volume densities nn >0.1 cm“^ 
(Figure!^, which are already ruled out by He Il/Lya. 
To conclude, in the context of our simple model, both 
high-ion metal-line absorption and He II and C IV emis¬ 
sion paint a consistent picture of the physical state of the 
gas. 

For completeness, we also searched for metal-line ab¬ 
sorption along the companion quasar ‘QSO b’ sight¬ 
line in our Keck/LRIS spectrum (resolution R ~ 1000 
and S/N ~ 60 pix“^). We do detect strong, satu¬ 
rated C IV absorption with Vciv > 10^'* "^ cm“^ and 
z = 2.2601. This implies, however, a velocity offset of 
« —1700 km s“^ with respect to the systemic redshift 
of the UM287 quasar, and thus from the extended Lya 
emission detected in the slit spectrum of Figure 2. Given 
this large kinematic displacement from the nebular Lya 
emission, we argue that this absorption is probably not 
associated with the UM287 nebulae, and is likely to be 
a narrow-associated absorption line system associated 
with the companion quasar. This is further supported 
by the strong detection of the rarely observed N V dou¬ 
blet. The large negative velocity offset —1370 km s“^ 
between the absorption and our best estimate for the 
redshift of QSOb z = 2.275 (from the Si IV emission 
line) suggests that this is outflowing gas, but given the 
large error ~ 800km s“^ on the latter, and the unknown 
distance of this absorbing gas along the line-of-sight, we 
do not speculate further on its nature. 

Finally, note that at the time of writing, there is no 
existing echelle spectrum of UM287 available, although 
given that this quasar is hyper-luminous r ~ 17, a high 
signal-to-noise ratio high resolution spectrum could be 
obtained in a modest integration. Such a spectrum would 
allow us to obtain much more sensitive constraints on the 
high-ion states G IV and N V, corresponding to Vciv < 
10 ^^ cm“^ and Vnv < 10 ^^'® cm“^, respectively, and 
additionally search for O VI absorption down to A^ovi < 
10^^ cm“^. If for example C IV were still not detected at 
these low column densities, this would raise our current 
constraint on nn by 0.5 dex to nn // 10cm“^ as shown 
in Figure Furthermore, the detection of metal-line 
absorption (at a velocity consistent with the nebular Lya 
emission) would determine the metallicity of the gas in 
the nebula, and Figure [^suggests we would be sensitive 
down to metallicities as low as Z ~ —3, i.e. as low as the 
background metallicity of the IGM (e.g. Schaye et al.j 
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5.2. Comparison to Absorption Line Studies 

It is interesting to compare the high volume densi¬ 
ties (nn > 3 cm“^) implied by our analysis to inde¬ 
pendent absorption line measurements of ga s densities 
in the CGM of ty pical quasars. For example Prochaska 
& Hennawi (20091 used the strength of the absorption in 
the collisionally excited C II* fine-structure line to ob¬ 
tain an estimate of riH — 1 cm“^ at an impact parameter 
of R± = 108 kpc from a foreground quasar, comparable 
to our lower limit obtained from the He II/Lya ratio. 
However, photoionization modeling of a large sample of 
absorbers in the quasar CGM seem to indicate that the 
typical gas densit ies are much lower ~ 0.01 <C 1 cm“^ 
( Lau et ^|2015 ), although with large uncertainties due 
to the unknown radiation field. If the typical quasar 
CGM has much lower values of nn 0. 01 1 cm~^ an d 

column densities of Ah ~ 10^° cm“^ (|Lau et al. 2015), 
this would explain why quasars only rarely exhibit bright 
Lya nebulae as in UM287. Indeed, eqn. § would then 
imply SBLyct = 5.4 x 10“^° ergs“^ cm“^ arcsec“^ in the 
optically thin regime, which is far below the sensitivity 
of any previo us searches for extended emission around 
quasars (e.g. |Hu &: Cowie 1987 Heckman et al. 


Christensen et al. 2006), altho ugh these low SB levels 


1991b 


may be reachable via sta cking (Steidel et al. 2011 Ar- 


rigoni Battaia et al.|2015 ). In this interpretation, quasars 


exhibiting bright SB ^ 10“^^ ergs“^ cm“^ arcsec“^ giant 
Lya nebulae represent the high end tail of the volume 
density distribution in the quasar CGM, a conclusion 
supported by the analysis of anot her giant nebula with 


discovered in the C 

luasars Probing Quasars survey 

Hen- 

nawi & Prochaska 

2013 

). In this system joint moc 

eling 


sightline piercing the nebular gas indicate that cool gas 
is distributed in clouds with i? ~ 40 pc, with densities 
riH — 2 cm“^, very similar to our findings for UM287. 

Absorption line studies of gas around normal galaxies 
also provides evidence for small-s cale structure in their 
circumgalactic media. Specifically, Crighton et al. (2015) 
conducted detailed photoionization modeling of absorb- 
ing gas in the CGM of a Lya emitter at z ~ 2.5, and 
deduced very small cloud sizes < 100 — 500 pc, although 
with much lower gas densities {n^ ~ 10“^ — 10“^ cm^) 
than we find around UM287. In addition, there are mul¬ 
tiple examples of absorption line systems at z ~ 2 — 3 
in the literature for which small sizes R ^ 10 — 100 pc 
have been deduced (iRau ch et ar]|1999[ |Simcoe et al. 
2006 Schaye et al. 2 0071, although th e absorbers may 
be larger at z ~ U.z (Werk et al. 2014). Also, compact 
structures with r ~ 50 pc have been directly resolved 
in high-velocity clouds in the CGM of the Milky Way 
(Ben Bekhti et al. 112009). Given their expected sizes and 


masses, such small structures are curre ntly unresolved 
in sim ulations (see discussion in § 5.3 of Crighton et al. 


6. MODEL SPECTRA VS CURRENT OBSERVATIONAL 
LIMITS 

In order to assess the feasibility of detecting other 
emission lines besides Lya from the UM287 nebula, and 
other similar extended Lya nebulae, e.g. around other 


quasars, HzRGS, or LABs, we construct model spec¬ 
tra using the output continuum and line emission data 
from Cloudy. In Figure P we show the predicted me¬ 
dian spectrum for the nebula at 160 kpc from UM287, 
resulting from our modeling. Specihcally, the solid black 
curve represents the median of all the models in our 
parameter grid which simultaneously satisfy the con¬ 
ditions 5.5 X 10“^® ergs“^ cm“^ arcsec”^ < SBLya < 
8.5 X 10“^® ergs“^ cm“^ arcsec“^, such that they pro¬ 
duce the right Lya emission level, as well as the emission 
line constraints He II/Lya < 0.18 and C IV/Lya < 0.16 
implied by our spectroscopic limits. Following our discus¬ 
sion in the Appendix, this grid also includes models with 
a harder (softer) aEuv = —1-1 (ctEUV = —2.3) quasar 
ionizing continuum, in addition to our hducial value of 
aEUV = —1.7. The gray shaded area indicates the max¬ 
imum and the minimum possible values for the selected 
models at each wavelength. 

For comparison we show our Keck/LRIS ia sensitivity 
limits from 1J2] calculated by averaging over a 1 arcsec^ 
aperture and over a 3000 km s“^ velocity interval (solid 
red line), together with the 3(7 sensitivity limits for 10 
hours of integrati on with the Multi Unit Spectroscopic 


Explorer (MUSE) (Bacon et al. 2010 solid blue line), and 


(Sharpies et al.| 

2006 

gold, orange, and dark-red solid 

lines), on the V 

L'i', computed tor the same spatial and 


spectral aperture. Note that these sensitivity limits can 
be lowered by assuming a certain amount of spatial aver¬ 
aging, following the relation S'Hiimit = SBi^ where 
A is the area in arcsec^ over which the data are averaged. 
Indeed, we employed this approach in ^ and averaged 
over an area of 20 arcsec^ to obtain a more sensitive con¬ 
straint on the He II/Lya and C IV/Lya line ratios, and 
this lower SB level is indicated by the red dashed line in 
Figure!^ In contrast with a longslit, integral-held units 
like MUSE and KMOS, as well as the upcoming Keck 
Cosmic Web Imager (KCWI, Morrissey et al.|2012 ), pro¬ 


vide near continuous spatial sampling over wide areas 
and are thus the ideal instruments for trying to detect 
extended line emission from the CGM. Thus for MUSE 
and KMOS, we have assumed that we can average over 
an area as large as 300 arcsec^, as shown by the colored 
dashed lines, and indeed this appro ach has already been 
used with the Cosmic Web Imag er (iMartin et a l. 2014 a I 
to study lower SB Lya emission ( Martin et al.||2014b ). 

Given these expected sensitivities, m I'lgurel^we indi- 
cate the principal emission lines that may be detectable 
(vertical green dashed lines), whose observation would 
provide additional constraints on the properties of the 
emitting gas. The large range of metallicities in our grid 
Z = 1O“®^0 to Zq^ results in a correspondingly large 
range of metal emission line strengths, whereas the Hy¬ 
drogen Balmer lines and He II are much less sensitive to 
metallicity and thus show very little variation across our 
model grid. 

Focusing first on the primordial elements, we see that 
He II is the strongest line, and in particular it is stronger 
than Ha. Indeed, if the Helium is completely doubly 
ionized then He II/Ha ~ 3, and although it decreases 
to lower values for lower ionization parameters (higher 
densities), it alwws remains higher than unity. As we 
have argued in ^ a detection of He II can be used to 
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Figure 9. Predicted median spectra for the models in our grid that satisfy simultaneously SBuya ~ 7 X 10 ergs ^ cm ^ arcsec 
He Il/Lya < 0.18, and C IV/Lyo < 0.16. The gray shaded area indicates the maximum and the minimum possible value for the selected 
models at each wavelength, showing the range of all the possible values, including the variation of the EUV slope, i.e. obuv = —2.3, —1.7, 
and —1.1 (see Appendix). Our Keck/LRIS 3cr sensitivity limit calculated in 1 arcsec^ and over 3000 km s“^, is plotted as a solid red line, 
together with the 3cr sensitivity of MUSE and KMOS (YJ, H, K gratings) for an exposure time of 10 hours (other colored solid lines). The 
red dashed line indicates our 3fT sensitivity limit average over an aperture of 20 arcsec^ (see ij^, while all the other dashed lines show the 
sensitivity averaged over an aperture of 300 arcsec^, i.e. SSiimit = SBi^/y/A. The principal emission lines are indicated by the green 
vertical dashed lines. The lines that may be detectable in the future, given appropriate physical conditions (i.e. nu,Z) in the targeted 
nebula are He II, [C III], C IV, [Si III], [O II], [O III], H/3, and Ha. 


measure the volume density nn of the emitting gas. Fur¬ 
ther, by comparing the morphology, size, and kinematics 
of the non-resonant extended He II emission to that of 
Lya, one can test whether resonant scattering of Lya 
plays an important ro le in the structure of the nebula 
( Prescott et al.|2015a| ). Naively, one might have thought 
that Ha would be ideal for this purpose given that it 
is the strongest Hydrogen recombination line after Lya. 
However, our models indicate that for photoionization by 
a hard source, the He II line is always stronger than Ha, 
and given that He II is in the optical whereas Ha is in 
the near-IR, it is also much easier to detect. 

Figure ^ shows that deep integrations in the near- 
IR with KMOS will consistently detect the Hydrogen 
Balmer lines Ha and H/3. When compared to the Lya 
emission, these lines would allows one to determi ne the 
extinction due to dust ( Osterbrock fc Ferland|2006 ). Fur¬ 
ther, at the low densities we consider (nn <C lO"* cm“^), 
any departure of the ratios Ha/H/3 and Lya/H/3 from 
their case B values provide information on the impor¬ 
tance of collisional excitation of Lya, which is exponen¬ 


tially sensit ive to the gas temperature (Ferland & Oster- 


brock 

19851. 

excita 

tion IS 


19851. In other words, the amount of collisional 


gas, which is set by the balance of heating and cooling. 
Photoionization by a hard source will result in a char¬ 
acteristic temperature and hence ratio of Lya/H/3 set 
by the ionizing continuum slope, whereas an additional 
source of heat, as has been postulated in gravitat ional 
cooling radiation sce narios for Lya nebulae (e.g. 
dahl fc Blaizot|[2012 1, would increase the amount o: 


Ros- 

__ col- 

lisionaily excited Lya and hence the ratio of Lya/H/3. 


Figure shows also that one could probably detect 
metal emission lines depending on the physical condi¬ 
tions in the gas, which are parameterized by nn and Z. 
In particular, if the gas has metallicity Z > Q.IZq, a 
deep integration with MUSE would detect C IV, [C III], 
and, for metallicity close to solar, also [Si III] A1883. In 
the near-IR, we see that a deep integration with KMOS 
would detect [O III] for Z > O.IZq, and [O II] for metal¬ 
licity close to solar. Note that for similar bright nebulae 
at different redshifts, it would be possible to detect other 
lines in extended emission for particular nn and Z com¬ 
binations, e.g. Si IV A1394, and [N IV] A1480. 

According to Figure [^ a good observational strategy 
is thus to look for the He II line, which appears to be 
the strongest and easiest line to detect, and our analy¬ 
sis in indicates that its detection constrains the gas 
properties to lie on a line in the nn-Z plane (see panel 
‘a’ in Figure [^. Following our discussion of C IV (panel 
‘b’ of Figure^, the detection of any metal line would 
define anotheiHine in the nw-Z plane, and the intersec¬ 
tion of these curves would determine the nn and Z of 
the gas. These conclusions will be somewhat sensitive to 
the assumed spectral slope in the UV (see Appendix), 
but given the different ionization thresholds to ionized 
Carbon to C IV (47.9eV), and Oxygen to O III (35.1eV) 
or O II (13.6eV), it is clear that detections or limits on 
multiple metal lines from high and low ionization states 
would also constrain the slope aEuv of the ionizing con¬ 
tinuum. 

To summarize, our photoionization modeling and anal¬ 
ysis provide a compelling motivation to find more bright 
nebulae by surveying large samples of quasars and 
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HzRGs, and conducting NB emission line surveys of 
LABs over large areas. Armed with the brightest and 
largest giant nebulae like UM287, one can conduct deep 
observations with IFUs, and combined with suitable spa¬ 
tial averaging, this will uncover a rich emission line spec¬ 
trum from the CGM and its interface with the IGM, 
which can be used to constrain the physical properties of 
the emitting gas, and shed light on physical mechanism 
powering giant nebulae. 


7. CAVEATS 

In section ^ under the assumption of photoioniza¬ 
tion by the central QSO, and in the context of a simple 
model for the gas distribution, we showed how our up¬ 
per limits on the He II/Lya and C IV/Lya ratios, can 
set constraints on the physical properties of the cool gas 
observed in emission. However, this simple modeling is 
just a zeroth-order approximation to a more complicated 
problem which is beyond the scope of the present work. 
In what follows we highlight some issues which should be 
examined further. 

Radial Dependence: for simplicity we have evaluated 
the ionizing flux at a single radial location for input into 
Gloudy. We have tested the impact of this assumption, 
by decreasing R from 160 kpc to 100 kpc, and find that 
our lower limit on the density increases by 0.4 dex. This 
results from the fact that the He II/Lya ratio varies with 
ionization parameter U, and our upper limit on the line 
ratio sets a particular value of U. By decreasing R, the 
density nn corresponding to this specific value of U thus 
increases by a factor R^. The variation of the ionizing 
flux with radius, should be taken into account in a more 
detailed calculation. 


Slope of the I onizing Continu um: we have assumed 
aEUV = —1.7 (Lusso et al. 2015). However, estimates 


for aF.TTv in the literature vary wide ly ( Zheng et al.|1997 
Scott et ar]|2004 Shull et ^|2012 1 , niost likely because 


of uncertainties introduced when correcting for absorp¬ 
tion due to the IGM or because of the heterogeneity of 
the samples considered. Furthermore, the shape of the 
ionizing continuum near the He II edge of 4 Rydberg is 
not well constrained. For detailed analysis on the sen¬ 
sitivity of our results to the ionizing continuum slope, 
see the Appendix, where we consider two different ion¬ 
izing slopes, i.e. Oeuv = —1.1 and —2.3. We find that 
a harder ionizing slope oeuv = —1.1 moves our lower 
limit on the density from nn ^ 3 cm“^ to uh ^ 1 cm“^. 
Thus, the uncertainty on the ionizing slope has an order 
unity impact on our constraints of the volume density. 
As discussed at the end of l|6l the detection of additional 


metal lines with a range ofionization thresholds would 
further constrain oeuv- 


Covering Factor: Based on the morphology of the 
emission we argued fc ^ 0.5, but assumed the value of 
/c = 1.0 for simplicity The fc drops out of the line ra¬ 
tios (see eqn. ^ and Jsl), however our model depends on 
fc, since we were selecting only models able to reproduce 
the observed Lya SB, which varies linearly with covering 
factor. We estimate that lowering the covering factor to 
fc = 0.4, only change our lower limit on the density at 
the 15% level. As discussed in section f|5 lowering fc 
results in a reduction of the number of models which are 
able to reproduce the observed Lya SB, because models 


with high nuNn valuse become optically thick, and thus 
over-estimate the Lya emission. In particular, there are 
no models which reproduce the observed Lya SB for low 
covering factors (/c < 0.3). Thus our conclusions are 
largely insensitive to the covering factor we assumed. 

Geometry: we have assumed the emitting clouds are 
spatially uniformly distributed throughout a spherical 
halo. This simple representation would need geomet¬ 
ric corrections to take into account more complicated 
gas distributions, such as variation of the covering factor 
with radius or filamentary structures. However, these 
corrections should be of order unity, and are thus likely 
sub-dominant compared to other effects. 


Single Uniform Cloud Population: our simple model 
assumes a single population of clouds which all have 
the same constant physical parameters Ah, u-h, and Z, 
following a uniform spatial distribution throughout the 
halo. In reality one expects a distribution of cloud prop¬ 
erties, and a radial dependence. Indeed, Binette et al.| 
(1996) argued that a single population of clouds is not 
able to simultaneously explain both the high and low 
ionization lines in the extended emission line regions of 
HzRGs, and instead invoked a mixed population of com¬ 
pletely ionized clouds and partially ionized clouds. While 
for the case of extended emission line regions (EELRs) 
around quasars, which are on smaller scale R < 50 kpc 
than studied here, detailed photoionization modeling of 
spectroscopic data has demonstrated that at least two 
density phases are likely required: a diffuse abundant 
cloud population with nn ~ I c ni~^. and much rarer 
dense clouds with rtn ~ 500 cm~^ dStockton et al.|[200J 


Fu & Stockton 2007| Hennawi et al.|^U9[ ). Further, these 


clouds may b e in pressure equilibrium with the io niz- 
ing radiation (Dopita et al. 2002 Stern et al. 2014), as 
has been invoked in modeling the narrow-line regions of 
AGN. Future detailed modeling of multiple emission lines 
from giant nebulae, an alogous to previous work on the 
smaller scale of EELRs ([Stockton et al.|2002||Fu & Stock- 
ton|2007[ ), might provide information on multiple density 


phases. 


In order to properly address the aforementioned issues, 
the ideal approach would be to conduct a full radiative 
transfer calculation on a three dimensional gas distribu¬ 
tion, possibl y taken from a cosmolog ical hydrodynamical 
simulation. Gantalupo et al. (2014) carried out exactly 
this kind of calculation treating both ionizing and res¬ 
onant radiative transfer, however this analysis was re¬ 
stricted only to the Lya line. Full radiative transfer cou¬ 
pled to detailed photoionization modeling as executed 
by Gloudy would clearly be too computationally chal¬ 
lenging. However it would be interesting to introduce 
the solutions of I-D Cloudy slab models into a realistic 
gas distribution drawn from a cosmological simulation. 
This would be relatively s traightforward for the case o f 
optically thin nebulae (e.g. van de Voort & Schaye 2013). 


8. SUMMARY AND CONCLUSIONS 

To study the kinematics of the extended Lya line and 
to search for extended He II A1640 and C IV A1549 emis¬ 
sion, we_obtained_deep_S£^troscopy of the UM287 neb¬ 


ula (Gantalupo et al. 


2014) with the Keck/LRIS spec- 


trograph. Uur spectrum of the nebula provides evi¬ 
dence for large motions suggested by the Lya line of 
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FWHMgauss 500 km s“^ which are spatially coher¬ 
ent on scales of ~150 kpc. There is no evidence for a 
“double-peaked” line along either of the slits, as might 
be expected in a scenario where resonant scattering de¬ 
termines the Lya kinematic structure. 

Although our observations achieve an unprecedented 
sensitivity in the He II and C IV line (SBao- ~ 
10 “^® ergs“^ cm“^ arcsec“^, average over 1" x 20"and 
Av = 3000 km s“^) for giant Lya nebulae, we do not 
detect extended emission in either line for both of our slit 
orientations. We constrain the He II/Lya and C IV/Lya 
ratios to be < 0.18 (Scr), and < 0.16 (3cr), respectively. 

To interpret these non-detections, we constructed mod¬ 
els of the emission line ratios, assuming photoionization 
by the central quasar and a simple spatial distribution of 
cool gas in the quasar halo. We find that: 


• if the gas clouds emitting Lya are optically thick 
to ionizing radiation, then the nebula would be 
~120x brighter than observed, unless we assume 
an unrealistically low covering factor, i.e. fc ^ 
0 .02, which is in conflict with the smooth mor¬ 
phology of the nebula. Thus we conclude that the 
covering factor of cool gas clouds in the nebula is 
high fc > 0.5, and that the gas in the nebula is 
highly ionized, resulting in gas clouds optically thin 
(A^hi < 17.2) to ionizing radiation. 

• The He II line is a recombination line and thus, 
once the density is fixed, its emission depends pri¬ 
marily on the fraction of Helium that is doubly ion¬ 
ized. On the other hand, the C IV emission line is 
an important coolant and is powered primarily by 
collisional excitation, and thus its emission depends 
on the amount of Carbon in the C’*'® ionic state. 
As we know the ionizing luminosity of the central 
quasar, and the Lya SB of the nebula, constraints 
on the He II/Lya and C IV/Lya ratios determine 
where the gas lives in the nn — Z diagram. 

• Photoionization from the central quasar is consis¬ 
tent with the Lya emission and the He II and C IV 
upper limits, provided that the gas distribution sat¬ 
isfies the following constraints: 

a) uh 3 cm“®, 

b) Ah 10 ^° cm“^, 

c) i? < 20 pc. 

If these properties hold through the entire nebula, 
it then follows that the total cool gas (T ^ 10"^ K) 
mass is Me < 6.4 x 10^° M©. 


Because the Lya surface brightness scales as SBLya oc 
tihAh, whereas the total cool gas mass as Me oc Ah, 
observations of Lya emission cannot independently de¬ 
termine the cool gas mass and nn (or the gas clu mping 


factor C), wh i ch lim ited the previous modeling by Can- 


talupo et al. (2014). Our non-detection of He Il/Lya 


combined with photoionization modeling allows us to 
break this degeneracy, and independently constrain both 
riH and Me- 

Our results point to the presence of a population of 
compact {R < 20 pc) cool gas clouds in the CGM at ISM- 
like densities of nn ^ 3 cm“® moving through the quasar 


halo at velocities ~ 500kms“®. It is well known that 
even by 2 ~ 2, the gas in the massive M ^ 10®^ ® Mq 
halos hosting quasars is expected to be dominated by 
a hot shock-heated plasma at the virial temperature 
T ~ 10^ K. Gool clouds moving rapidly through a hot 
plasma will be disrupted by hydrody namic instabilities 
on the cloud-crushing timescale (e.g. Jones et al. 1994[ 
Schaye et al. 1 120071 \ ” 

2015 


.ye et al.||2007i lAgertz et aLl|2T)07| [Crighton et al. 
f Scannapieco fc Bruggen||2015p 


tcc«l-3Myr 


— 

20pcy V 


20pc/ V500kms“i 


-1 


/ »cl/»halo V^^ 

^ 1000 J 


( 10 ) 


where we assume that the Lya line trace the kinemat¬ 
ics of the cool clouds, and that the cloud-halo density 


contrast is of the order of 1000 (rihaio ~ 10“® cm“®). 
If there is hot plasma present in the halo, these clouds 
are thus yery short liyed, and can only be transported 
^ 0.7 kpc before being disrupted. These yery short dis¬ 
ruption timescales thus require a mechanism that makes 
the clumps resistant to hydrodynamic instabilities, such 


as confinement by magnetic fields (e.g. McClure-Griffiths 
et al.||2010||McCourt et al.||2015|, otherwise the popula- 
tion of cool dense clouds must be constantly replenished. 
In the latter scenario, the short liyed clouds might be 
formed in situ, yia cooling and fragmentation instabili¬ 
ties. If the hot plasma pressure confines the clou ds, this 
might compresses them to high enough densities ( Fall &: 


Rees 1985 [Mailer fc Bullock|2004| Mc) fc Miralda-Fscude 


19961 to explain our results. Emission line nebulae from 
cool dense gas has also been obseryed at the centers of 


present-day cooling flow clusters (Heckman et al. 1989 


McDonald 


et al. 


2010 |), albeit on much smaller scales 
5Ukpc. Ihe giant Lya nebula in UM287 might be a 
manifestation of the same phenomenon, but with much 
larger sizes and luminosities, reflecting different physical 
conditions at high-redshift. Detailed study of the hydro¬ 
dynamics of cool dense gas clouds, with properties con¬ 
sistent with our cons traints, moying through hot pl asma 
are clearly required (j^annapieco & Briiggen||2015D. 

As we showed in ® deep obseryations ID nr) of 
UM287 and other giant ne bulae with the new integral 
field units such as M USE (Bacon et al. 2010), KCW I 


(Morrissey et al. 2012), and KlVlUb ( [Sharpies et al.|2006 ) 


combined with spatial ayeraging, will be able to detect 
extended emission from other lines besides Lya (see Fig¬ 
ure . In particular, the strongest line will be He II 
which should be routinely detectable, and following our 
analysis, will enable measurements of the yolume den¬ 
sity uh of the gas. Specifically, a 10 hour MUSE inte¬ 
gration would correspond to a sensitiyity in He II/Lya of 
^0.01 (3(7 in 300 arcsec^ ), which would allow us to probe 
gas densities as high as nn = 1000 cm“®. Although we 
have argued that the UM287 is powered by photoioniza¬ 
tion, which is compelling given the presence of a hyper- 
luminous quasar, a non-detection of He II in a lOhr 
MUSE integration would imply such extreme gas den¬ 
sities in the CGM, i.e. nn > 1000 cm“®, that one might 
need to reconsider other potential physical mechanisms 
for powering the Lya nebula which do not produce He II, 


such as cold-accretion (e.g., Haiman et al. 


2000 Furlan- 
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etto et al.| 

2005 

IDijkstra et al.||2006al IFaucher-Giguere 

et al.|2010i 

Rosd 

ahi & Hlaizot|2012|), star-formation (e.g.. 

Gen & Zhen^ 

2U13D, or superwinds (e.g., Taniguchi & 

SE^pUDO' 

ianiguchi et al.|2001 

Wilman et al. 

20 ^). 


matics of the nebula in He II and Lya can be used to 
test whether resonant scattering of Lya photons is im¬ 
portant. Although Ha could also be used to test the 
impact of resonant scattering, it is always fainter than 
He II and redshifted into the near-IR, where a detection 
of extended emission is much more challenging. 

In a photoionization scenario, a 10 hr observation of 
UM287 or a comparable nebula with MUSE (or KCWI) 
and KMOS would result in a rich emission line spectrum 
of the CGM, which, depending on the properties of the 
gas (i.e. nn and Z), could yield detections of Lya, [N IV], 
Si IV, [Ne IV], C IV, [C III], [Si III], [O III], [O II], H/3, 
and Ha. This would enable modeling of the CGM at a 
comparable level of detail as models of H II regions and 
the narrow and broad-line regions of AGN, resulting in 
comparably detailed constraints on the physical proper¬ 
ties of the gas. 

Current estimates suggest t hat ^ 10 — 20% of qu asars 


exhibit bright giant nebulae (Hennawi et al. 2015) like 
UM287, and our results provide a compelling motivation 
to expand current samples by surveying large numbers of 
quasars with instruments like MUSE and KCWI. At the 
same time, this same survey data would enable one to 
compute a stacked composite CGM spectrum of quasars 
which do not exhibit bright nebulae, constraining the gas 
properties around typical quasars. 
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APPENDIX 


EFFECTS OF A VARIATION OF THE EUV SLOPE OF THE INPUT SPECTRUM 

We test the robustness of our results to the change of the slope of the EUV as ment i oned in Section §4.3| In particular, 
we model the extremes of the range allowed by the recent estimates of Luss o et al. (2015), i.e. oeuv = —1-7± 0.6. To 
fulfill the Qfox requirement of Strateva et al. (2005) as explained in Section J 513 the value oeuv = —2.3 and -1.1 imply 
at higher energies (30Ryd < hv < 2keVj a slope a = —0.36 and -2.93, respectively. In our fiducial input spectrum 
(oeuv = —1-7), the photoionization rate at the Lyman limit is 


-9 


1 L 

r = -— 5 - / —^a^du = 6.7 X 10 “ s 


(Al) 


while at 4 Ryd, i.e. at the ionization energy of He II, the photoionization rate is r 4 Ryd 1.0 x 10“^^ s“^. By changing 
the slope in the extreme ultraviolet from ttEUV = —1-7, to -1.1 and to -2.3, we increase the photoionization rate by 
^15% and decrease it by ^13%, respectively. Instead, for the same change, the Td^yd is increased/decreased by a 
factor of 2.6, respectively. As it is clear from the small changes in T, the Hydrogen ionization state is not affected by 
the change in slope, and the models are always optically thin. Conversely, as expected, the changes in TdRyd affect 
He II and C IV. The general trend is that a softer slope, e.g. oeuv = —2.3, produces fewer He II ionizing photons, 
and thus at fixed density the He III fraction will be lower, resulting in lower He II recombination emission. This thus 
leads to a lower He II/Lya ratio. Similarly, a softer slope is less effective in ionizing Carbon. In particular, at fixed 
ionization parameter U, the amount of Carbon in the C’*'^ phase is lower for a softer slope. 

In Figure [Tol we compare our grids of models with different EUV slopes at two different metallicities, i.e. Z = Zq, 
and 0.01.^0, m the He II/Lya versus C IV/Lya plot. The dependencies outlined above, are better visible in the plot 
for solar metallicity (upper panel) because the Lya line is mainly produced by recombinations and its behavior is not 
influencing the general trends. From the figure, it is clear that a grid with a softer slope (see grid with Oeuv = —2.3) 
can reach lower He II/ Lya ratios because the fraction of doubly ionized Helium is lower at high densities. In the same 
upper panel of Figure it is also evident that the simulation grids for different UV slopes all asymptote to a fixed 
He II/Lya ratio when Helium is completely doubly ionized, which occurs at slightly different nn (or equivalently U) 
for each slope. Note that the value of the asymptotic He II/Lya ratio varies slightly with slope. Indeed, as mentioned 
in section ^ since this asymptotic value is proportional to the ratio of the recombination coefficients of He II and 
Lya, the value depends on temperature (eqn. @)- Higher temperatures, which arise for a harder slope, lead to a lower 
asymptotic He II/Lya ratio. _ 

In the bottom panel of Figure [TOl we show the same comparison at Z = O.OIZq. In this case the trends are masked 
by the Lya line, which is powered also by collisions. Indeed, the saturation in the He II/Lya ratio is not appreciable 
because, given the dependence on density of the collisional contribution to the Lya line, the ratio is progressively 
lowered at higher density. However, it is still appreciable that the C IV/Lya ratio is moved to lower ratios for higher 
slopes above logU ~ —1.5. This is mainly due to the fact that Carbon goes to higher ionization state, lowering the 
fraction of Carbon in the C"*"^ species. Thus, in our case study, where the input spectrum is not well known, the 
dependence of the amount of C"*"^ on the slope of the EUV makes the C IV line a weak metallicity indicator. 

Changes in the slope oeuv only slightly modifies the constrain ts o n tt-h that we previously obtained. In particular, 
since the He II/Lya ratio gives the stronger constraints, in Figurej^we show how a variation in the EUV slope affects 
the selection of uh (compare Figure 10 and [TT|). This Figure highlights in green the parameter space favored by our 
upper limits (the lines show the location of the upper limit He II/Lya = 0.18). The mild change in the location of the 
line is explained by the dependencies outlined above. At a fixed low metallicity, where the Lya line is an important 
coolant, i.e. logZ < —1.5^0, a harder slope moves the lower limit boundary implied by our measurement on the 
He II/Lya ratio to lower densities. Indeed, the expected increase of the He II line due to a harder slope is washed out 
by the increase in the emission in the Lya line due to collisions. Thus, our constraint on the density that we quote in 
the main text is weakened from uh ^ 3 cm“^ to rin ^ 1 cm“^. On the other hand, at higher metallicities, a harder 
UV slope will doubly ionize Helium at higher density, moving the lower limit boundary implied by our measurement 
to higher densities. For example, at solar metallicity, the limit is moved to > 100 cm“^ from > 40 cm“^. 

Thus, in conclusion, our ignorance on the slope of the EUV has a small effect on our density constraints and makes 
the C IV line a weak metallicity indicator. However, as discussed at the end of the detection of multiple metal lines 
with a range of ionization energies would indirectly constrain aEuV) and simultaneously constrain the metallicity of 
the gas. 
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Figure 10. Hell/Lya versus CIV/Lyo log-log plot. Our upper limits on the Hell/Lyo and CIV/Lya ratios are compared with the Cloudy 
photoionization models that reproduce the observed SEnyc* ~ 7 X 10“^® ergs”'^ cm“^ arcsec“^. Upper panel: comparison of the model 
grids for different EUV slopes (obuv = ~1-1) ”1-7, —2.3) at Z = Zq. A harder obuv completely doubly ionize Helium at higher density. 
Bottom panel: same as the upper panel, but at Z = O.OIZq. In this case, the Lyo line is also powered by collisions, reshaping the 
trajectories (see text for explanation on the trends in this Figure). In both panels, the models are color coded following the ionization 
parameter U, or equivalently the volume density nn (see color bar on the right). The green shaded area represents the region defined by 
the upper limits of the UM287 nebula. See Figure [ll] for a better visualization of the constraints on the physical parameters. 
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log^[Z0] 

Figure 11. Schematic representation on how a variation in qieuv affects the constraints in riH and Z. The green area highlights the 
region of the parameter space selected by the upper limit He Il/Lya < 0.18 (see panel ‘a’ of Figure [7||. The solid, dashed, and dotted lines 
show the location of this upper limit for qeuv = —1.1 , and —2.3, respectively. It is evident uhat a change in the ionizing slope do 

not affect our main conclusions. Namely, if the nebula is photoionized by the UM287 quasar, there should be a population of dense cool 
gas clumps with very small sizes (<tens of pc). 













